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ABSTRACT 



The basic Mitrovic's method for analysis and design of linear feed- 
back control system is studied. Some further properties of the basic 
Mitrovic's equations and curves are investigated and applied to analysis 
and design of linear feedback control systems. The method developed in 
this paper gives more flexibility to the basic Mitrovic's method. 
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1. Introduction. 

The Mitrovic's method is an analysis and design technique ©f linear 
feedback control systems using algebraic equations and graphical methods „ 
The method is based on conformal mapping of the S-plane into the real 
domain through the characteristic equation of a feedback control system s 
and the basic principles involved in the Mitrovic's method are derived 
by applying the basic theorem by Cauchy in the mapping procedure. 

The fundamental theories are adequately developed and applications 
of the method are demonstrated in references 1 and 2. The main advantages 
of the Mitrovic's method seem to lie in the fact that all operations are 
carried out in the real domain and that it has potentiality of controlling 
or adjusting multiple variable parameters. 

The analysis of a system in the Mitrovic°s method is done by analyz- 
ing the relative location of the Mitrovic's working point M with respect 
to Mitrovic's curves s and design of a system is to reshape the configura- 
tion of Mitrovic's curves and the M point associated with the system of 
interest so as to meet the design specifications. The analysis and 
design techniques are deduced from the fundamental properties of the 
Bo vs curves for 5 = 5l demoted by \~t >i with the M point as listed 
below. 

1) If the M point is located at o^so^t on the Bq vs Bi curve 
for £ = $i then the characteristic equation associated with 
the I7i curve and M point has a pair of complex conjugate 

roots T s -^COtu tJi-3* &ni 

2) If tangent lines are drawn from the M point to the Bq vs B^ curve 

for ^ = 1.0 (the ITo curve ) s then the negatives of the slopes of 
the lines are the real roots of the characteristic equation. 

1 



3) If the M point is located on a region enclosed by a si curve 
encircling in the clockwise direction, and if the M point is 
in the first quadrant on the Bq vs plane 3 then all of the 
roots of the characteristic equation have 5 gtaater than ^ 1 . 

In this paper further properties of Mitrovic's equations and curves 
on the Bq vs Bx plane are investigated to extend the method to problems 
which are not treated in the literature , and to give the method a better 
flexibility. 

In Section 2 the curves are studied to provide a procedure to 
evaluate frequency response of closed and open loop systems on the 
Bq vs Bx plane so that one can obtain information for frequency response 
and time response simultaneously. Constant band width curves are intro~ 
duced on the Bq vs B^ plane to design a system to meet the band width 
specification precisely with other specifications such as ^ and Kv 
specifications . 

In Section 3 fourth order system curves are studied to provide 
fourth order system charts to fit any fourth order system. The analysis 
and design procedures developed in references 1 and 2 can be worked out 
with preconstructed curves on the charts and also a different approach of 
design technique is tried with fourth order charts. 

In Section 4 an analytical design technique is derived so that the 
labor required in constructing Mitrovic's curves can be avoided or mini- 
mized. 

Thus each section can be regarded as independent from another; 
however, the principles developed in Section 2 may go with any of the 
sections to give information on frequency response. 
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2. Frequency responses of closed and open loop systems on the Mitro- 
vic's Bq vs B^ plane. 

Mitrovic's curves for 5 = the I ~o curves , are obtained by mapping 
the Imaginary axis of the S-plane in to the Mitrovic's plane. 

For a characteristic equation of a system 

Rs) — S’V&tt-i S 7 ' ■ - - + S + & 0 = o ( 2 -i> 

Mitrovic's B^ vs B^ curve for is defined by 



*1 , 4-1 

B =H 4* fart) co* 

K~2- 



where Cd n and £ are the natural frequency and damping rati© 
respectively of a root of the closed loop system,, and 
(finis') is a function of ^ only. (See Appendix I.) 
which is a curve on the Mitrovic's Bq vs plane independent of the last 

two coefficients &q and a^ of the characteristic Equation (2-1). 

Mitrovic's working point M on the Bq vs B| plane is defined by the 
coefficients a^ and a^ such that 

M-MCaj.a,) 

thus, for a given characteristic equation, a set of Bq vs curves and 
an M point are uniquely determined and conversely. 

For a particular value of damping coefficient ^ s ® the equation 
defining the Incurve is, (using proper values of ^(0) and replacing 
CO^byO) the frequency on imaginary axis of the S-plane in Equation (2-2)), 

B = a,co - irco 4 -t-a 7 Go 6 -- • • • 

“ - , (2-3) 

B 0 =a,a) a*£ o*+d*ar- * - • 
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Because of the nature of the lo curve. the frequency dependence of 
Equation (2-1) can be represented by the lengths of vectors from the H 
point to the to curve. This indicates that frequency response of a 
closed loop system can be evaluated on the Mitrovic”s plane. 

In this section relations between the frequency response (magnitude 
and phase) of a closed loop system and the Ho curve with the M point on 
the Mitrovic's Bq vs Bj_ plane are investigated. Band width of a closed 
loop system is evaluated and the locus of the M points for a system to 
meet the specified band width is constructed. 
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2-1. Interpretation of equations (magnitude of frequency response) on the 
Bq vs B l plane. 

Consider a closed loop transfer function of a feedback control 
system 



_C_ , _A/(s) 

K { F(S) 



A t(s ) 



(2-4) 



S'”-* 3*1 S w + ... +cL,S+a, 

where N(S) is a polynomial of S of order less than n. 

Note that F(S) r 0 is the characteristic equation from 
which the associated curve and M point are determined. 

The magnitude of the frequency response of the system represented by 
Equation (2-4) is 

I Nfi/tti i 

(2-5) 



|c , | [N(jq))| 



and I4H)“’} can be represented in terms of &q s Bq and 
a^, B^ respectively where B^ and Bq are the Mitrovic’s variables defining 
the r 0 curve which are functions of CO defined by Equation (2-3). 

Rewriting F(S) in Equation (2-4) as the sum of even F(S) and odd 



F(S), 

F(s> = (a 0 + • • ) 

then F( j 03 ) = Re~( F (j « 3 )} ■+- j I rr\\ F ( j 

=[d 0 -(a a c^-d. <f of +-■■)} + -( a 3 w *-& s U) 

therefore by using Equation (2-3) 



fe{Rjot>y - a 0 -(<3*c o*-a*co“+ • = d 0 —Bo 

IwjFCj'to')} = oo[3i-(a 3 af— 3 s co%- -)j = cu(<a,-Bj) 

thus, Equation (2-5) can be written as 




(2-6) 



5 



The interpretation of Equation (2-6) in vs plane is as fol- 
lows . The magnitude A at frequency CO is the ratio of the length of two 

vectors N(jO) ) and (a Q - B Q ) +j COCaj^ - B^). 

Since a^ and slq are the coordinates of the point M(a^ a^) on the 
B 0 vs plane, and Bj. and Bq at any frequency are determined by the co- 
ordinate of a point at o) on f7 » the numerical value of the denominator of 
Equation (2-6) can be determined easily. 

Suppose the 17 curve and M point configuration is as shown in Fig. 

2-1 and choose a point P on the 17 curve at which the frequency is oi , 
then the coordinates of P are (B^(«o), Bq(o))). Locate a point Q at 

(a^, Bq(co)) so that Q is at the same ordinate as P and is on the line 

B 1 s a l » then the distances from M to Q and Q to P are 

piQ-|a o -B b «»')| , Qtp = ja,-|3, Cal’) | respectively. 

Locate a point R on the line segment QP such that QR = 0)'<2.p and asso- 
ciate a frequency oi to the point R, then Q.R= o)*Qp= cO ja,-B|(0o)| 
and MR=(MO* ■+ )’ fe = | ( a 0 - +j oi (d,- B|(o)Y) J which is the 

value of the denominator of Equation (2-6) at frequency (0 . 

As 01 varies, the point P moves along the 17 curve and the point Q 
moves along the B^ = a^ line maintaining the same ordinate as P. The 
point R will describe a curve which is represented by the dotted curve 
and denoted by the OOCa^-B^) plot. The 0)^ - B^) plot is a well defined 
curve if the 17 curve is well defined and has the following common proper- 
ties for all systems: 

1) It has frequency on the curve as a parameter such that two 

points constituting a pair of points one on the 17 curve and the 
other on the C0(a| -B^> plot have the same frequencies with the 
same ordinates. 
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Fig. 2-1 Illustration of obtaining Frequency- 
response for a^syfetem with no- zeros. 
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2) It starts (or the originating point is) at (aj,0) and passes 
through the intersection point of the curve and the 

B^ s a^ line. 

3) The points for 0)= 1 on the 17 curve and Q)(a^ - B|) plot are 
identical. 

4) It lies between the 7 curve and the B^ = a^ line for CL)<1 
and dually. 

If any circle centered at the point M passes through a point at 
tO r on the £0(a^ - B^) plot with radius of r^ then r^ represent the 
value of the denominator in Equation (2-6) at ££>= The circle cen- 

tered at the M point is called the magnitude circle. The numerator of 
Equation (2-6) can not be evaluated from the I© curve as the denominator. 
It is possible to draw separate Mitrovic's curve for N(S) and obtain 
N(jcO) but this is labirious if N(S) is of high order. 

The simplest case is for N (jcO ) independent of frequency. This 
is the case for a system that has no zeros in its closed loop transfer 
function. 

Suppose a closed loop system has no zero, then 

A _ a ° 

|ra 0 -e | or 

|(a 0 -BJ+jcuf«5L r B l )| - 

The magnitude A^ at 00 = 0^ is a^/r^ and conversely if a magnitude 
A : is given then the frequency to give A : A^ is found by drawing a 

magnitude circle with radius equal to a Q/Ap and locating the point or 
points where the circle intersects with CD(a^ - B^) plot. 

Some particular points in the frequency response of a closed loop 
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system such as the resonant point, zero db point and half power point 
(band width) can be found very simply. 

Drawing a magnitude circle tangent to the u)(a^ - B^) plot gives 
the resonance point. Magnitude circles with radii of Sq and fl aQ give 
the 0 db point and band width respectively. 

For a system with N(S) in Equation (2-4) of order other than zero, 
the procedure is slightly modified but still the basic idea is the same. 

This is illustrated with Fig. 2-2. Modification must be made to 
take care of the frequency dependence of the numerator in Equation (2-6). 
Since the denominator in the equation defining the magnitude of fre- 
quency response bears the same form regardless of the order of N(S), the 
construction of the £t)(a^ - B^) plot is the same for all systems. 

Rewrite Equation (2-6) as 

(2 . 7) 

and interpret as follows; 

The quantity at any 0) on the left side of Equation' (2-7) is deter- 
mined by drawing a magnitude circle, and it is represented by the length 
of radius of the magnitude circle. On the right side of Equation (2-7), 

| f\|(j&>)| has some finite value for every a) . Then a number which 
multiplies to j N (j tb) j to give the equality in Equation (2-7) is the 
reciprocal of the magnitude. 

In Fig. 2-2, the value of |[S/(jCo)| is represented on a line drawn 
from M point in an arbitrary direction. This line is denoted by |N(j6o)| 
and at each point on the line, there associate cD such that if the point 
N is at Cx) on the line then 

MA/=|Nyw)| and A=FTF)/r 
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Fig. 2-2 Illustration for obtaining frequency response 
for a system with zeros. 
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' 








Conversely, o) = GO-^ for given magnitude A ; Aj is determined by search- 
ing such Cd that for two magnitude circles with radii of r and Ar cross 



quency. This can be done by a few tries and errors. 

The arguments developed above are not restricted by the order of 
system equation. The only restriction is to have curves plotted on 
the Bq vs plane with the same scales of B^ and Bq axes to provide 
simplest graphical solutions. Bq vs B^ curves can always be plotted on 
such axes scales by using suitable frequency scaling technique. 

As mentioned previously, the higher the order of the numerator 
polynomial of a system, the more complicated the method is. The problem 
is to provide a convenient method to represent the frequency dependence 
of the numerator of a closed loop system equation. 

In the following paragraphs, such methods are presented for the 
cases of first, second and third order numerators. 



the OMa^ - B^) plot and the 




line respectively at the same fre- 
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2-2. Representation of frequency dependence of N(S) cn the Bq vs Bj plane. 

In the previous section, it is seen that the frequency dependence 
of the numerator polynomial N(S) in a closed loop system equation must 
be represented as a distance from the M point. 

In Fig. 2-3 it is illustrated for N(S) of first order. When a 
closed loop system has one zero then N(S) is of the first order as 

N(S)=b,S 4-a 0 an d J SJ (j co) | — | d. 0 -+- j bjCi) | 

The real term of N(jCO ) is independent of frequency. By construct- 
ing a rectangle with a fixed side of length equal to a^ and another side 
which is in right angle with the fixed side varying linearly with C 0 , 
j f\} Cjo3) j can be represented by the length of the diagonal. 

In Fig. 2-3, a point T 2 is located on the Bq = a^ line such that 

mt 2 = b i* the li ne segment MT 2 is calibrated linearly with CO such 

that CO z 0 at M and CO = 1.0 at T 2 , and if the same calibration is ex- 
tended beyond the point T 2 , then at eaeh point on the line Bq = 3q to 
the right of M, there associate 00 linearly with the distance from M to 
each point. For any point at CD= Cdj., MT^ = b^ £0^ . Noting that 

TT^ = aQ, the length of diagonal of the rectangle MT^Ta^ is equal to 

MT = |a 0 +j | n ( j fOi) | 

For a system with two zeros in its closed loop system equation, N(S) is 
of second order polynomial as 

|s|(S)=b^+-bi s + ^0 and J N £/ co)j= b,C0 | 

The imaginary term can be represented in exactly the same way as in case 
of first order N(S). 

The real term is now also frequency dependent, and can be repre- 
sented on the B^ = a^ line with quadratic calibration. 
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In Fig. 2-4, bj4) is represented as in Fig. 2-3. To represent the 
2 

term (a^ - b^Od ) the line * a^ is calibrated as follows. The point 
PQ(a^,0) is labeled GJ = 0. From the point Pq up each point on the 

r a^ line is calibrated with such that if a point Pj on the s 

_ 2 

line reads 0) = then the distance = bjw^ . Noting that 

MP 0 = % ■ 

MfJ-FTPo-fFi- a.-M* 

By forming a rectangle MT^TP^ , the length of the diagonal MT is 

MT — | (a e - j b, oi | = | n Cj £o)| 

For a system with N(S) of third order or higher, |n(jCo)| can be repre- 
sented in a similar way as in case of second order N(S) with some modifi- 
cations as the case of first order N(S) modified to account the second 
order N(S); however, it requires laborious calculations. For example, 
for third order N(S) 

N(S) = b^S 3 b,s 4 a 0 (2-8) 

N(jcO ) = (a o -b 2 w z ) 4j(tjC0- b 3 C0 3 ) 

which indicates that the Bq = aQ line must be calibrated to account 
3 

(b^O) - b 3 Ca> ) which requires cubic calculations, and quadratic calcula- 
tions for fourth order N(S) and so on. 

Those calculations can be avoided by applying Mitrovic's basic 
equations for the curve to Equation (2-8) and by proceeding the same 
way as was done to represent F( jcd ) by means of CO (a^ - B^) plot. 

This is illustrated with Fig. 2-5. 

Since N(S) is of third order, the 1^ curve obtained from Equation 
(2-8) is a straight line and this is denoted by . Let Mitrovic's 
variables associated with Equation (2-8) be Bm and Bq^. Then 
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f^u is defined by 

B on -^ 



(2-9) 

The point Mjj(b^ , Sq) which is analogous to the point M(a^ s Sq) for 
the characteristic equation F(S) r 0, is determined by the coefficients 
b^ and aQ of Equation (2-8). 

If E^ n and B qn axes are chosen identically with B^ and Bq axes* then 
has a slope of b 2 /b^ and is located on a horizontal line passing 
through M on the Bq vs B^ plane. 

If 0) (b^ - B Ln ) plot is constructed from l^ N and the line B^ s B^ N = 
b^ , then by the same reasoning as for G)(a^ - B^) plot, for any fre- 
quency &) L the quantity | W (j 0\) | = \(a o - E^ + j ©j ( ^ - B JN ) | 
where B^ N and B QN are coordinates of a point at 6)s CO^ on , is 
equal to the distance from Mjj to T 2 where the point T 2 is at 6)= on 
the 4)^ - B 1n ) plot. 

Since | N(jC0 )| and | (3q - Bq) -t-J CO (a^ - ) j must be compared to 

get A the magnitude, and since N(jo3 ) is measured from M N whereas the 
latter is measured from M, it is necessary either the 0)(b^ - E^ N ) plot 
is translated horizontally to the right or the C0(a^ - Bj_) plot to the 
left so that Mjj and M coincide. In Fig. 2-5, the 6)(b^ - ®ijj) plot is 
translated to the right by amount of a^ - b^ , then 

|n(j'«i)|= MJi =Tfr 3 

The method representing |n(JC 0 )| for N(S) of up to third order is 
presented. For high order N(S), |n(j&> )| can be represented by an analo- 
gous method as in case of third order N(S); however, it is inevitable to 
be laborious to represent N(j‘C0 ) as the order of N(S) increases. 
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On the Bq vs plane for feedback control systems with adequate 
relative stability , the M point is usually located in the region near to 
the origin as compared with whole region covered by the Fo curve. This 
is due to the fact that curves for higher ^ are biased towards the 
origin. Therefore in obtaining frequency response usually a first part 
of the curve is required. 

It must be noted that, since the basic Mitrovie's equations are not 
restricted to the characteristic equation ©f a feedback control system 
but apply to any polynomial with positive coefficients, the principles 
developed above apply to any transfer function which is expressed in the 
form of a quotient of two polynomials provided the Bq vs curve for 
s 0 and the M point associated with the transfer function are defined, 
Therefore it is obvious that the magnitude in the frequency response of 
an open loop system can be evaluated in a similar way on the same 
Bq vs B^ plane where the magnitude of the closed loop system is evalu- 
ated. 

In general the open loop frequency response can be obtained by con- 
structing the rj curve and by locating the M point associated with the 
open loop transfer function and then by proceeding exactly the same way 
as for the closed loop system. In Section 2-7, some examples for evalu- 
ating the phase of the frequency response for both open and closed loop 
systems are shown, which will well illustrate the above arguments. 
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2-3. Determination of band width of a closed loop system on the 

Bg vs plane. 

For any given system, band width defined at magnitude of 1//5 or 
half power point can be determined precisely by the method described in 
Section 2-2. 

For a system with no zeros, band width is found at the intersection 
point of CO (a^ - B^) plot and a magnitude circle with radius equal to 
/2a 0 . 

For a system with zeros, band width can be found by trying a set of 
2 concentric magnitude circles with radii of r and fz r such that the 
two circles cross the CO(a^ - B^) plot and |n(JcO)| respectively at a 
same value of frequency. 

Frequently in the study of feedback control system the point of band 
width alone of the frequency response is of interest and also sometimes 
upper bound and/or lower bound of band width is of interest. 

By the property of the equation defining band width and associated 
geometry on the Bg vs B^ plane it is possible to set up upper and/or 
lower bound of band width for a given system. Among the various points 
in a frequency response curve of a given system,, if band width only or 
approximate value of band width is required, then by setting up the lower 
and/or upper bound of band width, the effort in plotting £>(a^ - Bj_) can 
be economized or even without plotting any portion of 6d(a^ - R^) the 
approximate value of band width can be determined. 

In this section the three cases for a system with (1) no zeros, 

(2) with one zero, (3) with two and three zeros are considered. 

For a system represented by Equation (2-4), and using Equation (2-6), 
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the band width defined at half power point is the frequency to satisfy 
the equation 

|C fnj-J — |HQ°»I _L 

+jo) (<Si-3|)| 

or I ( a o - I3 0 ') +jw(a, -B|')| -fz I N(jc*))| ( 2 - 10 ) 

All that has to be done is to interpret Equation (2-10) on the Bq vs B| 
plane. The left side of Equation (2-10) can be interpreted by the 
60 (a^ = Bj) plot as shown in Section 2-1, and for the right side is 
needed a modification to the method described in Section 2-2 to account 
for the factor [2. 

Consider a closed loop transfer function with no zero. 

Let 6)^ be the band width, then from Equation (2-10) 

| +jco t> (a | -B J (o ) b )|=/2 ( 2 -ii) 

must hold, where B^( C0 b ) and Bq( 6) b ) are the coordinates of a point on 
the V 0 curve at C0= 60 b . 

Equation (2-11) implies that a magnitude circle with radius of 
J~ 2 aQ intersect with the C0(a^ - B]_) plot at C0 b which is described in 
Section 2-2. Since the V 0 curve tends to infinity in some direction as 
60 tends to infinity, ^(a^ - B^> plot also tends to infinity. Noting the 
finite radius of a magnitude circle and that the 0)(a^ - Bj_) plot starts 
from a point inside the magnitude circle with radius /T a^ it is seen 
that there exists at least one intersection point between the circle and 
the 60 (a^ - B^) plot. Depending upon the relative location of the M 
point with respect to the Vq curve and also depending upon the geometric 
shape of the fp curve there may be more than one such intersection point 
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implying more than one band width to exist. 

On the Bq vs plane with the curve, by drawing a magnitude 
circle of radius (~2 a q and by inspecting the values of 0) at all inter- 
section points of the circle with the curve if there are any 3 and by 
estimating the shape of the 60 (a^ - B^) plot 9 information about the band 
width can be obtained. Also the point at CO s 1 on and the points 
where the curve intersects with s a^ s Bq s Eq and Bq s (1 + 
lines provide additional information. 

Suppose a system is represented by the 1"^ curve and M point as in 
Fig. 2-6. A magnitude circle C Ab with radius of fl aQ is drawn. The 
points P^, P2* P3, P4, and P5 represent the points where the curve 
intersects with C Ab , Bq = a Q , Bq = (1 ■+ /~2)a^ or = a^ . 

Suppose at each point of those intersections * frequencies are read 
to be 00^, GO^, Cd^t Cd A and as shown. 

The latter part of the Yq curve is neglected since it is assumed 
that for that part of , the C0(a^ - B^) plot stays outside of C Ab 
In an actual situation this can be judged by looking at the quantity 
C0 b (a^ " b i(^ 6)) and 6)7 (a^ _ B (6)7)) where <5J b and 6)^ are the fre- 
quencies at the points where the latter part of Q intersects with 
Bq z ( 1 + /~2) a o and Bq = aQ respectively. 

The - B^) plot resulting from the portion of the 1"^ curve 

above the line Bq = ( l+/2)aQ can never cross C Ab , therefore 9 then 
there are four possibilities for band width as 



(1) 


C0 3 - 


— 60^ 


(3) 


(x> l — 63(, — 


(2) 


a** 


^3 


(4) 


£0 b ^ &>! 
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Fig. 2-6 Determining approximate band width 
for a system with no zeros. 1 
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Each of the four possibilities above can be evaluated by looking at the 
point at CO: 1 on the curve. 

Suppose &>3>1 then the <5)(a^ - B^) plot for 60 — CO3 stays outside 
of the 17 curve in turn outside of C Ab s indicating 60 b < 6)3 . 

Suppose W 3 <1 then Q3P3 >^(3^ - B ( 0) 3)) and the GO (a^ - B^) 
plot must cross C Ab somewhere between £0^ and £0^ therefore 
^<60 b <c0« < 1 OT 0) 3 < C0 b < 1 

Suppose fi): 1 is located between ^ aB <* ^3 or <co 0 s 

then, since the curves for £4)3 < £0< Co ^ are inside of C Ab * the 
60 (a^ - B^) plot must cross C Ab between 4) r 1 and 6)3 3 therefore 
1< 00 b < 0)3 • In this case there is still a possibility that a) b < 61)3 
but this can be checked by looking at the quantity 4)(a^ - B^) for some 
00 less than CO2 . 

Suppose 60^ < 1 < 6D2 » then the a) (a^ - B^) plot stays outside of 
C Ab for 1^60<C02 and 60^4)3 . The portion of the 4)(aj_ - B^) plot for 
6O2<<W<CO3 may be inside of C A 2 . If this is the case then &)2^^b < ^3 , 
again the possibility that £0 b < CO2 must be checked. If COCa^ - B^) for 
60 2 <: t0 < a)3 has no point inside C Ab then COCaj. - B^) for 60£ 1 is out- 
side C Ab and ^<1, If the points and P£ are very close, then 

1 = <W, = 0J Z 

Suppose 60^>1, then there is no possibility that 60^ is in the In- 
tervals ( 1 , 60 2 ) and 60> 6O3 . Since the 60(3^ - B^) plot at 60 = 0 is 

at the point (a^, 0) and the 60(a^ - B^) plot at 60= 1 is outside of 
C Ab , band width is definitely 0< 6D b < J. < COj and possibly 
00 2 < £ Oh <. 0O3. 

If either one of the frequencies 60^ and 6O3 is equal to 1, then 
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r 1 and other possibilities can be investigated as above. 

Thus the four possibilities listed above are analyzed. The geometry 
involved in Fig. 2-6 is so simple that the relative locations of 
63( a l - plot for 60*0^ , 6) 3 4 c0^g) 4 , 0) 3 , 

and 0* with respect to are determined almost by inspection 

using the point C0= 1 as a guide, resulting in determination of upper 
and/or lower bound of band width. 

It must be noted that, depending upon the system of interest, there 
will result various configurations of the curve and the working point 
M (a^,aQ); however, the basic principle of reasoning is the same as 
above . 

For a system with one zero, the band width defining Equation (2-10) 
is to have frequency dependent terms in both sides of the equation and 
requires some modification in the method applied for a system with no 
zeros. If a system has one zero in its closed loop transfer function, 
then the equation defining band width is 

|(a t> -B < ,(co^Vj«b( a .-B l (cOb^|=/2 la.+jMbl ( 2 - 12 ) 

where b^ is the coefficient of the term in the numerator of 
the given closed loop transfer function. b q( 
are the same as in Equation (2-11). 

Suppose the 1"^ curve and M point configuration is as shown in Fig. 
2 - 7 . To account for the frequency dependent term f2 b^c3 on the right 
side of Equation ( 2 - 12 ) the line Bq : ( 1 +J~1 ) aQ (or Bq = line) 
is calibrated with CO linearly such that the points Q and T^ correspond 
to 6)= 0 and 63=1 respectively, and the length QT]_ is equal to 
f~ 2 b^ . Then the distance from Q to any point at Gl)=G)j on the line 
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Fig. 2-7 Determination of band width for a system 
with one zero. 
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Bq = ( l-+f2 )ag will represent fl bj^j and it is obvious that 

MW s fZ |d 6 +jb,C0 , | which is the quantity on the right side of Equation 

(2-12) for (D r 60’ 

Suppose the C0(a^ - B^) plot is as shown by the dotted curve in 
Fig. 2-7. Noting (She nature of the cO(a^ - Bj_) plot that it always 
starts at the point (a^ 9 0), crosses at 0) s 1 and at the point 
where intersects with - a^ line* it is easy to guess or sketch the 

(DCa^ - B^) plot. 

The band width C0 b is some 60 such that Equation (2-12) holds. 

This can be interpreted in Fig. 2-7 as follows? Suppose a magnitude 

circle is drawn such that it pass through the point W at £Q = a>’ 

on the B 6 = (1 +V2) a 0 line and suppose it reads 60 " 

on the o^a^ - B^) plot. If CO* = uf then Equation (2-12) holds at 

frequency equal to CO , therefore CO^ = CO . 

To determine band width one must find a C A such that it crosses 

both &)(a^ - B^) plot and Bq - ( 1 •+ fl )a^ at the same frequency. 

For a C ^ if 0>V co then the two frequencies are either cO>Ci) 

or 6o’< 60 ” . By drawing another C Ac0 ’ until the inequality between 

CO and tO is reversed from the previous set of CO and CO the bounds 

of band width can be set up. 

For example, for the configuration of curves as shown in Fig. 2-7 
1 is assumed. 

For the two magnitude circles C a q and 

Frequency on Bo = (1 + /~2)a 0 line Frequency on 60 (ap - Bi ) plot 

C A0 60'r 0 < CO* » 0) 0 

C A i 60 = 1 > 60" = 4), 

in which the inequality sign is reversed,, therefore 0 < c0 0 < C0 fc < < I 
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It is obvious in Fig. 2-7 that there is no possible band width for 
63 < &)q . There may be another band width for 03 > 4)^ . This can be in- 
vestigated by trying some magnitude circles. 

C aQ in Fig. 2-7 is identical to C Ab in Fig. 2-6. If > 0 in Fig. 
2-7, then stays on C^q for all 00 , which is the case for a system 

with no zeros and agrees with Fig. 2-6. It is readily seen from Fig. 

2-7 that the effect of zero and forward gain of a system on band width 
can be analyzed. 

For two systems, one for no zero and the other for one zero, if the 
characteristic equations are the same, then band width of the first system 
is less than that of the second. This is confirmed in Fig. 2-7 as for 
the first system » 00 q and for the second system C0 b >(0 q . 

The effect of varying forward gain and the zero on band width for a 
system with one zero is seen from Fig. 2-7 by writing the numerator of the 
closed loop transfer function as 

N(S) s b L S+a Q = KS +• KZ 

Increasing either K or Z is to expand both C a q and which is to 

increase COq and fi3^ resulting in increased band width. Since K and Z 
do not enter into any coefficient other than a^ or of the character- 
istic equation, the ^ curve remains unaffected by either K or Z, 

Suppose K is fixed and Z is variable. Since a^ and b^ are independ- 
ent of Z, the C0(aj - B^) plot and the lines Bj^ = aj. 9 4 /2 b^ 

remain unchanged. Varying Z moves the M point and the line Bq « slq ver- 
tically, thus increasing Z expands C a q and C A ^ resulting in increased 
band width. 
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Suppose Z is fixed and K is variable, since a^ „ and bj. are 
linearly varying with K» all of the lines in Fig. 2-7 and the oKa^- B^) 
plot are moved. The M point is moved along a line with some positive 
slope and C a q and C A i are expanded with increasing K resulting in in- 
creased bandwidth. 

In a system which has positive coefficients in both numerator and 
denominator polynomials of its open loop transfer function*, b| is always 
less than a^ assuming unity feedback. This means that the bandwidth of 
a system with one zero is always less than that of a system which results 
from replacing b^ with a^. Therefore for a system whose ^ and M point 
configuration as shown in Fig. 2-7 bandwidth must be 
& 0 < < 0d 2 

where 0)2 is the frequency on the 60(a| - plot at which the £i)(a^ - B^) 
plot and magnitude circle C A £ intersect. C A 2 drawn such that it passes 
through T 2 where T 2 is located on the Bq : ( 1 + f~2 )a^ line at the dis- 
tance of {H a^ measured from Q. 

Above inequality is useful to set up the bounds of £0^ when is a 
variable and when it is required to locate M point to meet some required 
bandwidth specification. It must be noted that the bounds determined by 
the above inequality become poorer as the difference between bj and a^ 
gets larger. 

It is interesting to note that if b^ = then the open loop 

transfer function is of the form, assuming unity feedback, 

3|S+a 0 

^ •+ <2^., s'*" 1 + - • - +a.j 

which is a type 2 system, and this system has greater bandwidth than a 



27 



system whose open loop transfer function is 



r/c\ = 

-- - + a 2 s i -f(a,-b,)s 

For a system with two zeros 9 bandwidth can be evaluated in a simi- 
lar way as illustrated in Fig. 2-7. One modification is required in 

Fig. 2-7 to account for the frequency dependence of f2 (aQ- b2£0^) 3 

2 

where b£ is the coefficient of the S term in the numerator of the 

closed loop transfer function. This is done by calibrating the line 

% 

as in Fig. 2-4 to represent (aQ - b 2 &)^ ) by the distance from M 
to the point at (j) z on the = a^ line. T© account for the factor 
J~2 j, the origin of calibration or the 60s 0 point is chosen at a point 
on the B^ = a^ line below the axis s© that it is apart from the M 

point by a distance of f2 a q. Then the line is calibrated with 

Z 

60 such that any point at 60 = 60^ on the line will measure /2 
by the distance from the 60= 0 point (originating point) to the point 
where 00= 60^ . 

For a system with 3 zeros, the quantity J ( Hq - b2tO^)+j'ft)(bi»b260^) | 
which is the right side of the bandwidth defining equation can be repre- 
sented in a similar way as in Fig. 2-5. The only difference is to ex- 
pand the translated C0(b^ - B^ plot radially from the M point by factor 
of fl . 

Then for both cases of the above s bandwidth can be evaluated by 
applying the same analysis as for the system with on© zero in Fig. 2-7. 
The following numerical examples illustrate the arguments developed so 
far. 
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Example 2.1 

Obtain frequency response (magnitude) of the closed loop system 



JLfS'N- - T. 



(2-13) 



The system represented by Equation (2-13) is a third order system. 
Since it has no zeros the bandwidth will be less than the frequency at 
the point where the [~ Q curve and Bq = (1 +■ /~2 )a^ line intersect. If 



f”o curve needed is below the line Bq = (1 + / 2 )aQ = 9 . 65 . 

In Fig. 2-8 the frequency response of the system ( 2 - 13 ) is repre- 
sented on the Bq vs B^ plane. 

From the characteristic equation of the system (2-13) the \ 0 curve 



is defined by 

Bn=3of 



(2-14) 



which is a straight line with slope of 3 passing through the origin. 
The magnitude of the frequency response is defined by 

«■ 



A(io)= 






(2-15) 



where B and B ft are defined by Equation (2-14). 



“1 ” “0 ' 

The point M (2, 4) is located, and the lines B^ = a^ = 2, Bq = aQ = 4 
and Bq = (1 -+-/~2 )Sq = 9.65 are drawn as shown. The d0 ~ B ^) P^ ot is 

constructed as illustrated in Section 2-1. 

The <i)(a^ - B^) plot starts from the point (2, 0) and crosses the 
PJ curve at 0) r 1 and the point (2, 6) at which the ^ curve and the 
line B^ = 2 intersect. 

For 0 < 0) < 1 and for to >1.42 the 0)(a^ - B^) plot lies to the 
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right of the curve and for 1.0 < 10< 1.42 to the left. The point at 
00 s .5 on the CO(a^ - B^) plot is located on a horizontal line passing 
through the point at £0 = .5 on the curve and the distance from the 
line B| s 2 to the point at 0) : .5 on the &)(a| » B^) plot is ©ne-hal£ 
of the distance from the line B^ : 2 to the point at CO* .5 on the 
curve so that the distance from the ; 2 line to the point at 60s .5 
on the 4)(aj. - Bj.) plot is equal to 0) (2 - B^(C0 )) j ^ ^ which is the 

value of imaginary term at Od * .5 in the denominator of the equation 
defining the magnitude of the frequency response. 

Some points on the 60 (ai - B^) plot for some pertinent values of 
frequencies are determined likewise and those points are connected by a 
smooth curve to give the 60 (a^ - B^) plot as shown. Thus each pair of 
two points on a horizontal line, one on the curve and the other on the 
60 (aj. - Bj_) plot has one value of frequency and the distance from the 
point M (2, 4) to any point at frequency 0) on the £0(aj_ - B^) plot is 
the value of the denominator of Equation (2-15). 

Magnitude circles are constructed as shown. Noting that the numer- 
ator of Equation (2-15) is constant at ag r 4„ it is seen that the ratio 
of a^ s 4 to radius r of a magnitude circle which passes through a point 
at frequency of 0) on the fl)(a^ - B^) plot is the magnitude A 0 ) ) . 

At 60 = 0, the radius of magnitude circle is 4 which gives A = 1. 

As the frequency increases from zero, the radius of the magnitude circle 
decreases, in turn increasing A. At Cl) ~ 1.18 magnitude circle is 
tangent to the o)(aj_ - Bj_) plot and has radius r = . 75 which implies 
resonance peak Mp^ = 4/. 75 = 5.33 
at resonance frequency cO r = 1.18 
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For frequencies from CO = 1.18 on, the radii of magnitude circles in- 
crease, in turn decreasing A. 

A magnitude circle with radius of fl - 5.65 crosses the 
“ B^) plot at Ci) r 1.76 which implies the bandwidth is 
CD b = 1.76 

The frequency response up to Cl) = 1.8 is thus represented on the Bq vs 
B^ plane. 
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Example 2.2 



Given an open loop transfer function of a unity feedback servo as 

r , -v t 2 / Q&O (£ + 37 ) 

(X ^ S) S^ + A I)(g + lt>70) 

obtain the magnitude vs frequency plot of the frequency response. 

The closed loop transfer function is 



c,. J2.jo(>bg+wsooo 

s 3 +/ £jijg2+ J2267/&+ W8040 (2-16) 

Since Equation (2-16) deals with large numbers , it is not convenient to 
work on the Bq vs plane associated with Equation (2-16). 

By using a transformation or frequency scaling 
S = 100s 



Equation (2-16) transforms to 

_C v_ /g./s -f ___ 

I? S 3 + \e. r 7fs 2 ' + /2.2,67s + 

Then the equation defining frequency response is 

V I V V-3 +j/ 2 .i cc> t I 



(2-17) 



(2-18) 



where B^ t and B^ t are Mitrovic's variables of B^ vs B^ curve 
for $ = 0 associated with the system (2-17) and is re- 
lated with CO, the frequency associated with the system (2-16) 
by 

01 = 100 co t 

In Fig. 2-9 the V 0 curve, (^(a^ - B^ t ) plot and magnitude circles 
representing the numerator of Equation (2-18) are constructed. 

The frequency dependence of the imaginary term in the numerator of 
Equation (2-18) is represented on the B Qt - 2a Qt line to the left of 
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B^ t = a^ t line. (Note that the frequency dependence of the imaginary 
term in the numerator of Equation (2-18) is represented on a line dif- 
ferent from that illustrated in Section 2-1.) 

As illustrated in Section 2-2, the distance from the point 
M (12.267, 4.48) to a point at frequency 6) ± on the = 2ag t - 8.96 

line is equal to | 4 . 48 +j 12.267 C0 t | 

The magnitude circles are drawn so that at any frequency the numer- 
ator and denominator of Equation (2-18) can be compared. 

Since the M point is located at the point where the coordinates are 
not convenient numbers, reading off the radii of magnitude circles is 
not convenient. 

This can be simplified by providing a scale which is calibrated 
with the same scale of coordinate axes so that the radius of any magni- 
tude circle is read off easily. This is shown in Fig. 2-9, for example, 
at 6) t = .5, the value of |n( jcO t )| | ^ _ 5 is read to be 7 » 5 on the 
scale provided on the line B^ t = 12.267 by drawing a magnitude circle 
passing through the point at C0 t = .5 on the line Bg t = 2aQ t and simi- 
larly for the denominator | (a 0t - B 0t ) +j 6) fc (a lt - B 1{ .) j j^_ 5 s 6.0, thus 

the following table is constructed: 



35 




36 




' 










Example 2.3 

Suppose for the system in Example 2.2 only the bandwidth is of in- 
terest. The equ at ion 

*jcd ± (/a. 27-13, Jc0^)\=/z\wg+J/2./a) t \ (2-19) 

defines the bandwidth £0^ of the system in Equation (2-17). 

By drawing a magnitude circle with radius equal to 1. 414(4. 48+j 12. 1) 
which is the quantity on the right side of the equation above at = 1, 
and noting the frequencies in the neighborhood of the point on the rj, 
curve where the circle intersects with the V 0 curve, and also noting the 
relative location of the “ B^ t ) plot with respect to the V 0 

curve, some information can be obtained about the bandwidth. This is 
illustrated with Fig. 2-11. 

In Fig. 2-11, the curve is drawn for the system (2-17), the work- 
ing point M (12.27, 4.48) is located, and B^ t = a^ fc r 12.27 is laid off 
the same way as in Fig. 2-9. The line Sgt = (1 t 1.414)a 0t = 10.8 is 
drawn. 

From the intersection point of Bg t r 10.8 and B^ t - 12.27 lines to 
the left, a distance of 1.414 b^. r 17.1 Is measured off on the Bg t = 

10.8 line and labeled Cd t = 1.0. The radius of a magnitude circle pass- 
ing through the U = 1.0 point on the Bg t = 10.8 line is equal to 
| 1.414(4.48 -+j 12.1) | = |l.414 N(j 1) |. This circle denoted by C Al crosses 
the Yl curve between the points at &) t = 1.0 and 6) t = 1.1. 

It is easy to predict that 60 < 1 is not possible. If it is sus- 
pected that 6). < 1 is impossible, it can be checked as follows ; 
bt 

The check is based on the fact that, for the bandwidth, the magni- 
tude circle must cross the G0 t (a^ t " B Q t ) Pl ot an< * 1-414 b^ t scale 
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represented on the line Bg t = (l-t/2 )ag t = 10.8 at the same frequency, 

and the fact that if any magnitude circle is drawn and if it crosses the 

^t^ a lt " B 0t^ P lot and l*414b lt P\scale (B Qt . - 10.8 line) at 60 t = 0) t " 

and 0) t = 60 j. respectively with either 60<6o" or a/> a) , then for 

existence of a bandwidth frequency, there must be another magnitude circle 

to give a new set of CO and c0 with inversed inequality from the previous 

set of CO and 60 . 

The magnitude circle C , gives 
A1 

60” > i and d0 t ’ = 1 then 60 t < u)^. 

Since C A ^ is drawn with oj^= 1 and it crosses the curve between 
(x) t = 1.0 and 60 t = 1.1 and 1.0 <Q\< 1. 1. This is reasoned without 
drawing any portion of the ^t^ a lt " B lt^ P^ ot * 

Draw a magnitude circle with 60= 0. This circle is of radius equal 
to [2 ag t tangent to the line Bg t = (1 + fl )aQ t z 10.8 and represents 
f2~ I N<j 6> )| = [~2 a Qt by its radius. Since 60 t (a 1( . - B lt ) plot ap- 
proaches to the point at Ci) t = 1 on the curve from the point (12.27., 0), 

the circle must cross the k) t (a^ t - B^ t ) plot at some frequency greater 

than zero to give an inequality 

=0 y 60 1 >0 then c0 t ’ < 0)^ 

The inequality sign is not reversed from the previous set of and 
and this implies there might be no frequency such that, for 0 < I0 t < 1 , 
a magnitude circle crosses the 60 fc (a^ t - B^ t ) plot and the {~2 b^ t (i) t 
scale at the same frequencies; in other words, z 0)^ r 6i) bt ma y not 
be possible for the frequency range from zero to 1.0. However, there 

may be some frequency between 0 and 1.0 at which the inequality is re- 

versed since the interval of 60 1 from 0 to 1.0 is equivalent to a rela- 
tively large interval of original system frequency (0= 100^, 
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Locate the point at ^ r .5 on the " B^) plot which is 

half way from the line B lt = a lt = 12.27 to the point at 0> t s .5 on the 
[g! curve. This point is located inside the magnitude circle correspond- 
ing to 6i^= 0, therefore, no bandwidth is possible for (0 t ^* 5. 

For the magnitude circles corresponding to „7 < 1, the inequal- 

ity oS<(0 is not changed. This can be checked by comparing the dis- 
tances from the M point to the fg curve and the [2 b^j. ft) t scale and 
noting that the “ B^ t ) Pl°t H®s to the right of the ^ curve for 

.7 < a ) t < i.o. 

For example, a magnitude circle for r . 7 will cross the ^ curve 
at a point where is greater than .8 to give Ifi” >.8 and 0) t < o)J! 

Thus for the bandwidth, . 7 <£. 0i), < 1.0 is not possible. 

Dt 

The only part left to check is the interval .5 <{Jt) t < .7. Knowing 

that i , <> 

for CO t z' .5 Cx> ^ < CO and 

for (0’s .7 a)'<af 

it is quite sure that for the interval .5 < £0 t < .7 the same inequality 
will be maintained. If a magnitude circle for * .6 which is the mid 
point of the interval of the frequency from .5 to .7 is drawn, this cir- 
cle will pass the V Q curve at the point where (0 1 = .6 implying the 
point at = .6 on the ^ t ( a i t “ B^ t ) plot will be inside the circle 
in turn 60^.<G0 t which is the same inequality as the cases of 6t.' t = .5 and 
CO^ = .7. Thus the same inequality is maintained for &) t from 0 to 1.0 and 
1 is not possible. 

The tests made above are carried out without plotting any portion of 
the 60t( a lt " B lt) pl Qt except the point at Ci> t = .5 on the B]_ t ) 

plot is located. The tests can be carried out with a sketch of the 
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60 c ( a i t ~ B it) pl Qt and the sketch can be made simply with reasonable 
accuracy by noting that the horizontal distance from the line B^ t = a^ t 
r 12.27 to the (jO t (a^ t - B^ t ) plot is a fraction of that to the curve 
and the fraction is the value of 0 ) fc . 

For an interval of frequencies (&)^, 6 ) 2 ), if the inequality 
60 ’ < CO” (err a)' > a)") 

holds at both ends and then the inequality is maintained through- 
out the interval unless the ^ curve or the ~ B^ t ) P^ ot i s oddly 

shaped in the interval, or in other words, unless the rates of change of 
the distance from the M point to the - B^ t ) plot is not fluctuat- 

ing in the interval. 

In this example, if the C0 t (a^ t “ B ^ t ) P^ ot sketched, then one 
can see that the rate of change of distance from the M point to the 
^t^ a lt ~ B lt) Pl° c i s steadily increasing except in the region of small 
frequencies. 

Knowing that 61 )^ >1.0, the next thing to do is to find another 
magnitude circle to insure an upper bound of 61 )^ as illustrated in Sec- 
tion 2-3. 

Noting that crosses the fo curve at the point where the fre- 
quency is slightly less than 1.1 and also noting that the (i) t (a^ t ~ B^ t ) 
plot lies to the left of the 17 curve for ^>1, it is predicted that 
(0b t will not be very far from Q) t z 1.1. Therefore it is natural to 
try a magnitude circle which passes through Cl) t = 1.1 or = 1.2 on 
the /2 b^ t Ct) t scale. The magnitude circle which passes the point at 
£ 0 t = 1.2 on the line Bg t s . (1 *[2 )ag t is tried and denoted by C^2* 

The circle represents the right side of Equation (2-19) at 0) t - 1.2. 
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It passes the curve at a point where the frequency is less than 1.2 
which implies it will pass a point on the “ B lt) Pl Qt where fre- 

quency is less than 1.2, thus C0 t = 1.2 clearly indicates an upper bound 
of bandwidth and 

1.0 < Cd < 1.2 

b t 

To determine 6) bt precisely, the ^t^ a lt “ B lt^ Pl ot f° r l*0<A) t < 
1.2 is needed. 

A magnitude circle C^ b which passes through the point at Ol ) t = 1.1 
on the /~2 b lt tO t scale crosses the ^ t (a^ t - B^ t ) plot at cO t = 1.1 and 
bandwidth is determined to be 

This agrees with Fig. 2-10. 

Since Fig. 2-10 or 2.11 is for the system represented by Equation 
(2-17) which is related to the original system by S = 100s , the original 
system bandwidth (0^ is 

(0 b = 100 L 0 bt = 110 

By substituting s = j C0 bt = 1.1 and S = j 6i) b - HO into Equation 

(2-17) and (2-16) respectively, the magnitudes calculate to 




l c M 



I 



43 



2-4. Locus of M points on the Bg vs plane for constant bandwidth. 

In Section 2-3, bandwidth of a given system is determined. For a 

given system with fixed M point which is determined by the coefficients 

a^ and ag, bandwidth is determined by varying CO until some frequency 

(or frequencies) CO = cj is found to satisfy Equation (2-10). 
b 

Suppose, conversely, CO in Equation (2-10) is fixed at bandwidth 
6t) b , then there must be some a^ and ag to satisfy Equation (2-10), and 
the pair (a^ , ag) defines the M point on the Bg vs Bj plane for a system 
to have bandwidth equal to 4)^. If all of such pairs (a^ , ag) satisfying 
Equation (2-10) are plotted on the Bg vs B^ plane, then a curve which 
represents the locus of M points on the Bg vs B^ plane for fixed band- 
width Ci) b is obtained. This curve is called the constant bandwidth curve 
for bandwidth equal to ti)^ or merely constant £0 b curve. 

If a constant curve is constructed on the Bg vs B^ plane, and if 
the M point is chosen at any point on the curve, then the resulting sys- 
tem will have bandwidth equal to C0 b * 

Since the synthesis of a feedback control system by Mitrovic’s 
method is to choose a suitable M point to guarantee the required system 
performance specifications, the constant £0 b curve will provide informa- 
tion in choosing the M point to meet specified bandwidth requirement. 

By fixing 0J = in Equation (2-10), the constant (0^ curve is 

defined by 

| =/2 | N (j CO b )| (2-20) 

where B^ (&) b ) and Bg ( 60 b ) are the coordinates of the 
points on curve at it) = ft) b . 

On the left side of Equation (2-20), B^ ( <X> b ) and Bg(Cd b ) are fixed and 
those values can be read off on the curve. 
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On the right side of Equation (2-20), noting that N(S) is the numer- 
ator polynomial of the closed loop transfer function, |n( J^ b)! contains 
ag which is treated as a variable on the constant 6) b curve. Moreover, 
if N(S) is of order other than zero, it may contain which is the co- 
efficient of the S* term in N(S), and b^ is not independent of a^ which 
is the other variable on the constant curve. Thus, even though 
Equation (2-20) is a perfectly general expression for defining constant 
curves, it does not clearly indicate the geometrical shape of the 
constant curve on the Bq vs plane. 

In the following, the cases of zero order, first order and second 
order of N(S) are considered. 

N(S) of zero order . 

When a closed loop transfer function has no zeros, N(S) = ag and 
Equation (2-20) becomes 

| ( a 0 -B.K)) +j °> b ( a,- B, <w>) I = n <2. 0 

by removing the absolute sign, and manipulating 

(a„+B 0 (io k jf- <(VB,W) = * <2-2i) 

Since 6> b , B 0 (<i) b ), B p(^ b ) are fixed. Equation (2-21) indicates a hyper- 
bola with following geometrical properties on the Bg vs B^ plane; 
the center of symmetry at ( Bj (a)p) , 
the vertices at (B i (fd b ) J (~l±J5) 
major axis is j3j = Bj 

The two branches of hyperbola (2-21) lie symmetrically with respect to 
the Bg = -Bg(6L) b ) line. The upper branch lies above the B^ axis with 
minimum point at (B^((0 b ), .414 Bo(COb)) which is one of the vertices. 
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For a stable system, only the portion of constant £*) b curve which lies 
in the stable region is of interest. 

It is interesting to note that the locus of the vertex of the branch 
of the hyperbola which lie in the first quadrant is similar to the fj 
curve, and the locus of the center of the hyperbola is the mirror image 
of the C curve. 

Those can be seen as follows. As shown previously the vertex of the 
upper branch of the hyperbola is at (B^(Ci> b ), .414 Bq( C 0 b )) and the center 
point is at (B^O^), -B Q ( 4) b )). 

Considering C0 b as a continuous positive variable, then the locus 
of the vertex point is defined by 

X = B ± (sO 

Y = .414 BjCtO - ) 

and in like manner the locus of the center point is defined by 
X = B X C 0 l> ) 

Y = -B 0 (CO) 

where X and Y axes are identical to B^ and Bq axes respectively. 

The two sets of equations can be compared with the curve defined by 
B t = B t ( CO ) 

Bq = B 0 (CO) 

Consider the consistency of principles Involved in constant 60^ curve 
and in determining bandwidth described in Section 2-3. 

In Fig. 2-12, suppose a constant C0 b curve is plotted from Equation 
(2-21) assuming the f? curve as given. The point P is at Cl) z CO. on 

O D 

the ^ curve. The constant C0 b curve is denoted by C The vertical 

axis denoted by Y is on the B^ = B^(cO b ) line which passes through the 
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Fig. 2-12 Showing consistency of Section 2-3 and 
Section 2-4 about bandwidth. 
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point P and V is the vertex. The horizontal axis and the lower branch of 
the hyperbola which lie below the B^ axis are not shown. 

The portion of the constant 6) b curve which is of interest for a 
stable system is the portion of C^ b under the curve. 

Suppose an M point is chosen arbitrarily on C ^ b as shown. Then a 
(a^ - B^) plot results as described in the previous sections and a 
magnitude circle C A with radius equal to (l a^ must pass through the 
point Q at which 0)- 6) b on the ^ (a^ - B^) plot to have a consistency 
for the principles involved in the constant 4) b curve and determining 
bandwidth as described in Section 2-3. This is to show that 

mq ./a a, 

Locating a point R at the intersection point of the B^ = a^ line and 
B 0 = B 0 <0) b ) line, 

1 / 2 . 

Ra«jfiR% r* [f p, (c*$ 

But (a^ , ag) satisfies Equation (2-21) and it rearranges to 

and HQ. -/I a 0 

N (S) of first order . 

When a closed loop system has one zero 
N(S) = b,SH* 

For unity feedback control system, b^ enters into a ^ and aj and b^ are 
not independent of each other. 

Suppose an open loop transfer function of a unity feedback control 
system as 

^ ^ b £> + 21 o 

"♦--•td.s 
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where and a^ are variable and d^S are fixed constants. 

c, . b,g +a 0 

Then -^{S'i s , ^ + . . . + (d| + b ,i S + a 0 

c /r , N b,£ + 

or R (s)- s *> + ^ + a, s t a, (2 ' 22) 

where = d n _j, , ... a 2 = ^2 and a l =d l' hb l 

N(S) can be written in the form of 

N(S) = (a L - d L ) S+a Q 

With this N(S), Equation (2-20) is manipulated to 

(A+5>W\) + [a,+(B,(a) b V^di'|] (2-23) 

where D (a5 b , d , ) - 2 [ P* (0) b ) + cOf ( B, (o^ - d, f j + d?4. 



Equation (2-23) defines the constant curve for any system whose closed 
loop transfer function has one zero. 

Noting that for D(&> b ,d^)>0, Equation (2-23) is an ellipse centered 
at (2d^ - B^(fit) b ) , -Bq( &) b ) ) with semi-axes of lengths equal to 



VDr^b’di) and y'ixcO^^j) . (Horizontal and vertical axes respectively. 

For each d-. , as CO, varies there results a family of ellipses, 
i b 

It is seen in Section 2-3 that for a system represented by Equation 
(2-22), larger b^ results in larger bandwidth. 

For any unity feedback system assuming all positive coefficients in 
its open loop transfer function, b^ is always less than a^ and its upper 
limit is a ^ which is the case of a type 2 system, whereas its lower limit 
is zero which is the case of a system with no zeros. 

For the upper limit of b^ = a^ corresponds to d^ = 0, and letting 
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= 0 in Equation (2-23), the constant curve is then 



(a 0 + b„ (<oo f ■ + w l (■ a i + 8, W)f = ^ OV) +< 3 W) 



Equation (2-24) can be applied for a type 2 system with one zero in its 
closed loop transfer function. For any system with b^< a^ or d^> 0 , 
if the constant 0) b curve is plotted by Equation (2-24), and if the M 
point is located on the curve, then the system will have bandwidth less 
than Ct)^, thus Equation (2-24) defines a locus of M points on the 
Bq vs B^ plane for some upper bound of bandwidth for a system whose 
closed loop transfer function is represented in the form of Equation 
( 2 - 22 ). 

This will be illustrated with a specific system 



and this curve is plotted in Fig. 2-13 as shown. 

If the system had no zero, then constant 6 ^ curves would be defined 
by Equation (2-21) with B^ and Bq defined as above. For this case, con- 
stant curves are constructed for 6 ^ : .1 to 6 ) b = .5 as shown in 
Fig. 2-13. If M point is chosen at any point on C 5 which is a constant 
curve for 6 ^ b r .5 then the system 



c^ ( 0 a * —21 



(2-25) 



with associated constant £t) b curves as shown in Fig. 2-13. 
For the system (2-25) the 17 curve is defined by 





(2-26) 



will have bandwidth equal to .5. 
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The dotted curve which is a parabola is the locus of the vertices 



of the constant G3 b curves for the system (2-26). 

A constant £) b curve of G) b = 1.0 for a system 




(2-27) 



which is obtained by letting b^ = a^ in Equation (2-25) is plotted by 
Equation (2-24) and denoted by Q in Fig. 2-13. If the axes scales 
are made the same, the q would be an arc of a circle, since, for 
G) b = 1, Bq(1) z 0, B^( 1) = 1 then Equation (2-24) is 



which is a circle centered at (-1, 0) with radius equal to fl . 

If the M point is chosen at any point on q then the system 
(2-27) will have bandwidth equal to 1.0. 

Suppose the point M(.407, .089) is chosen for both systems (2-26) 
and (2-27). Since the point M (.407, .089) is on both C 5 and q , 
the systems (2-26) and (2-27) have bandwidth 6) b = .5 and = 1.0 re- 
spectively. 

For system (2-26) 



I C r A 


.OP? 1 


1 H 0* 5 )|- 


+ .0S9 



For system (2-27) 



Then the point M (.407, .089) for the system (2-25) guarantees that 
0.5 < eo b < 1.0 
for 0 < b^ < a^ 



a£ + (a l + i 
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N(S) of second order and up . 



When a closed loop system has 2 or more zeros 
N(S) z •+ b i S? ■+ S' ■+•'** 

If the coefficients b^S for k^2 are fixed (in order for a sys- 
tem to be analyzed or synthesized by Mitrovic’s method on the Bg vs 
plane, the coefficients a^S for k>2 must be determined and likewise 
for b^S ), then by proceeding in the same way as the case of N(S) = b^S 
-t- aQ, the equation defining the constant , curve, which is similar to 
Equation (2-23) with due modification to account for the coefficients b£, 
b^, . . . will be obtained and likewise for the equation defining constant 
C0 b curve of certain upper bound for bandwidth. 

It must be noted that for a high-velocity constant system b^ is very 
close to a^ since the velocity constant Ky is related to a-^, aQ and b^ by 



Therefore, for a high Ky system, the constant C0 b curve obtained by let- 
ting b^ s a^ will be a good approximation to the actual bandwidth. 

Also in lag compensation, a^ = b^ and as shown in references 1 
and 2, lag compensation problems are solved by letting a^ = b^ , conse- 
quently, the constant h) b curve obtained from a^ s b^ will be good in 
lag compensation problems. 
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2-5. Bandwidth of third order systems. 

It is shown in references 1 and 2 that the equations of third order 
systems can be normalized to give a universal Mitrovic's chart for third 
order systems, and this chart is constructed in references 1 and 2. 

In the normalized chart, is a straight line extending from the 
origin to infinity with slope of unity on the Bg n vs B^ n plane, and the 
frequency CO fcn in the normalized system is related to the frequency CO in 
the unnormalized system by 



00 = k 00 tT) 

where ^ is a constant, the normalizing factor, and sub- 
script 7) denote normalized quantity. 

Assume a third" of der system with no zeros in its closed loop trans- 
fer function: 

C_fr\ 

s 3 * + 

It normalizes to, by transformation S = a 2 S 



C 

j 3 + a 17 ,d+ a 07l 



(2-28) 



where d l71 = cL , / 3^ 



= 3.0 /4 



(2-29) 



Equation (2-29) defines the relation between coefficients of the 
normalized and unnormalized systems, consequently the Mitrovic's variables 
are related by Equation (2-29) between normalized and unnormalized sys- 
tems , 



namely B |71 =B,/aJ 



(2-30) 



54 



Applying Equation (2-21) to the system represented by Equation 
(2-28), the constant (j) b curve is defined by 

(<2-077 +3 (rn ( 6 W > ) + ^bn (^iyi — ZB 07) (tipnl (2-31) 

where is the bandwidth of normalized system and is re- 

lated to the unnormalized bandwidth C0 b by o) b s ^^bn' 

B^ n and BQ n are Mitrovic's variables for the ^7 curve in the 
normalized third order chart, and are related with unnormalized 
variables B^ and Bq by Equation (2-30). 

Equation (2-31) defines a family of constant curves on the 

B 0n vs B^ n plane if a set of 4> bn is supplied. 

If this family of constant G) bn curves are superposed on a third 
order chart, then it will provide additional information about the third 
order system. 

Fig. 2-14 shows the family of constant bandwidth curves for third 
order systems whose closed loop transfer function is represented by Equa- 
tion (2-28) on the normalized Bq vs B^ plane. Note that the portions of 
constant bandwidth curves defined by Equation (2-31) which lie in the 
unstable region are omitted. 

Fig. 2-14, the constant bandwidth chart for third order systems, can 
be superposed on Fig. 2-15 which is the third order Mitrovic's chart con- 
structed in references 1 and 2 to provide additional information concern- 
ing system bandwidth in both analysis and synthesis of third order sys- 
tems . 

To determine system bandwidth is to locate the M point of given 
system on the constant bandwidth chart and consult with the constant 
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Pig 2-15 3rd order chart (Normalized B 



vs B curves) 5" 





bandwidth curves. Also it is an easy matter to choose an M point such 
that it guarantees a given bandwidth specification. 

From Fig. 2-14, the effect of moving the M point (in the stable 
region) on system bandwidth can be analyzed as follows. 

Suppose a point M (ap n 9 a Q n ) is located in the stable region. For 
fixed a^, increasing results in increasing the bandwidth as is 
easily seen from Fig. 2-14. This agrees with common experience. 

For fixed a^, the effect of changing a^ n on bandwidth is dependent 
upon the value of a^. This can be analyzed by noting the slope of con- 
stant (jL) bn curves. 

It is seen in Section 2-4 that the locus of the vertex point of the 

family of constant bandwidth curves, the hyperbolas s is 
X ^ Bjn(a)t'n') 

Y = 0.414 B 071 (aV) 

where X and Y axes are identical to the and axes respectively 
and B ln <"tn> and BQ n ( i0 tn ) are the equations which define the (7 curve. 
Noting that the 17 curve for the normalized system is a straight line 
emanating from the origin with unity slope in the first quadrant, then 
the locus of the vertex points of the constant curves is a straight 

line with slope of .414 defined by 

B()n = “ 414B ln < 2 ‘32> 

on the BQ n vs B^ plane. 

Since a vertex point of a hyperbola is the point at which the slope 
of the hyperbola changes its sign, every constant curve lying in 

the stable region is such that the portion of the curve which is above 
(to the left of) the line defined by Equation (2-32) has negative slope 
and dually. 
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Dividing the stable region into two regions by the line defined by 



Equation (2-32), it is easy to see that if a^ is fixed s then by increas- 



region above the line and dually for the region below the line. 

If a closed loop system has a zero in its closed loop transfer 
function, then the normalized system equation is 

bmxJ -+ S-07? 



/S 3 +ss 3 + a. lv 4 + a. 0 r, 



(2-33) 



where b^ n is normalized b-^ , the coefficient of S* term in the 
numerator polynomial of the unnormalized system, and the rela- 
tion between b^ and b^ n is 



hm = bi/aj 

The arguments developed in Section 2-4 for a system whose closed loop 
transfer function has one zero applies to the third order system repre- 
sented by Equation (2-33). 

By letting b^ n = a ln in Equation (2-33), then the resulting constant 
bandwidth curves which are defined by Equation (2-24) in Section 2-4 
(with due account of normalized third order ) ares 



0i7i (2-34) 

where B. and B rt are Mitrovic's variables for the 0 curve 
In On 

on the normalized third order charts. 

The portions of the ellipses lying in the stable region are plotted 
in Fig. 2-16 which can also be superposed on Fig. 2-15 to provide infor- 
mation about bandwidth for any third order system with one zero in its 
closed loop system transfer function. 
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As illustrated in Section 2-4, Fig. 2-16 constitutes a bandwidth 
chart for a system represented by 

r . , + < 









2.,-n y6 -f 3 4? 



and a bandwidth chart for the upper bound of bandwidth for a system rep- 
resented by 

c (s) ^ * a ° 

K l ‘ ' ^• J + a i ^ = + a,s + a. 

or equivalently for the normalized system in Equation (2-33). 

The effect of moving the M point in the stable region on bandwidth 
can be analyzed in Fig. 2-16 analogously as in Fig. 2-14. 

Note that the ellipse defined by Equation (2-34) has a center point 
moving along the line BQ n = B^ n in the third quadrant as <D bn varies s 
and also note that the ellipse has negative slope in the first quadrant. 



is increased. 

In Fig. 2-17, 2-18, 2-20 and 2-21, the constant bandwidth charts. 
Fig. 2-14 and 2-16 are plotted on expanded scales. The third order chart 
Fig. 2-15 is correspondingly expanded in Fig. 2-19 and 2-22 which go with 
Fig. 2-17, 2-18 and 2-20, 2-21 respectively. 

The following examples illustrate determining bandwidths of third 
order closed loop systems by using the bandwidths charts developed in 
this section, and also illustrate a design of system to meet exact band- 
width by using the third order system bandwidth chart. 
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Pig. 2-19 3rd prder chart (Normalize^ B- vs B curves) 
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1 Fig. 2-21. Constant bandwidth curves 

(Normalized 3rd order system, b. 
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Example 2.4 

For the system in Example 2.1 

£ 3 + 3 + V 

the bandwidth is found to be Cd b = 1.76. Now one can find the bandwidth 
on the third order bandwidth chart. 

The given system normalizes to, by using the transformation S = 3s s 

iM = '±tl 

R ^ -t- ^ + .Z22^J 

and the normalized bandwidth 4) bn is related with the unnormalized band- 
width CO by Ci) = 3 4) . Since the system has no zeros, one enters into 

b b bn 

the chart in Fig. 2-14 or the expanded chart in Fig. 2-17 or Fig. 2-20. 
The M point of the normalized system seems the chart in Fig. 2-17 is 
suitable. By locating the point M(.222, .148) on the chart in Fig. 2-17, 
6Jb n is found to be .59 by interpolating 60 bn = .55 and c0 bn = .6 curves. 
Then the bandwidth of the original system is 

^h=303 t/rt =( 3 ')(.??)=/-'77 

which agrees with the result in Example 2.1. 

For the system in Example 2.2 

c J2IO0DS +4<S-<?0DO 

,2267 /S+ V-V-8000 

the bandwidth is found to be C0 b r 110 in Example 2.3. 

This system is third order with one zero in its closed loop transfer 

function. If the bandwidth is determined from the third order bandwidth 

chart for a, = b, then the bandwidth so determined must be not smaller 
In In 

than 110. Noting that a-^ = 122671 and bj^ = 121000 are not greatly dif- 
ferent, the bandwidth determined from the chart for b^ n = a^ n will be 
very close to the actual bandwidth 110. 
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The system normalizes to, through the transformation S = 1671.1s, 

jC, .V „ +.POO 9 & 

A^ ■+ .<$ 2 ■+■ . oU 3 3 A + , 000 f b 

Entering into the bandwidth chart in Fig. 2-21 with the normalized 
M(.0438, .00096) point, the normalized bandwidth is found to be 
03 ^ 7 ) = . o 7 

by interpolating between the two curves for bn = *05 and bn = .075. 
Then the unnormalized bandwidth is 

<o b =a= r co Ml = ^7/.0 6»7) = //7 

which is very close to the actual bandwidth £0 b = 110 and agrees with 
the arguments developed in Sections 2-3, 2-4 and 2-6. 
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Example 2.5 



Given a system block diagram below 



Jo 



S(S+0(£ + <!) 



+ IftS 



adjust K,K a and K t with the following performance specifications: 

1) 5 : .3 

2) K v = 10 

3) 0) b = 20 

For the system without feeding back derivative signals K must be 
set to 90 in order to meet the K v specification, and the system is at 
stability limit which can be seen from the third order chart. 

The closed loop transfer function with the derivative signals feeded 
back is 

5 

R S 3 -»-(l 0-hK Ka^‘ Z -+-(^ + K 

and normalizes to, by using the transformation S = 33 s = (10+-KK a )s s 

C-ox- 

^ } 4 3 + / ) ' 3 ‘ +- + Zen 

where; a ln - ( 9 + K K*) / 0 0 + K Ka.'f 

a 0 n = F/Oo+KK^i 

and a 2 r 10 +- KK a = normalizing factor. 

The velocity constant and bandwidth of the unnorraalized and normal- 
ized systems are related by 
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k v = (io^ Khra') (C V7I 
CO b =($ O + 

K v /<D b = K vri / ^br, = '/*- 

where the subscript "n" denotes the normalized system quantity. 
Thus the ratio of the velocity constant and the bandwidth C0 bn of 
the normalized system must be maintained at 1/2. 

In Fig. 2-22a the fo and G curves are duplicated from the third 
order chart (Fig. 2-15) and then some constant bandwidth curves are 
superposed from the third order system bandwidth chart (Fig. 2-14). 

In order for the system to have "$= «3 the M point must be located 
below the r.z curve in the first quadrant on the BQ n vs Bj n plane. For 
the normalized system K vn = a 0n/ a ln indicates that for an M point located 
anywhere on the normalized Bg n vs B^ n plane if a line is drawn from the 
origin to the M point, then the slope of the line is equal to the normal- 
ized velocity constant K^. Therefore for this problem the M point must 
be at the intersection point of a constant fr) bn = ^bn curve and a line 
drawn from the origin such that the slope of the line is (l/2)dt) bn and 
moreover that M point must be such that it insures S = .3. For any 
chosen constant &) bn = £0 b ^ curve the corresponding straight line with 
slope 1/2 Cl)h n can be drawn and then the line almost always intersects 
with the constant &) bn = ^ b n curve, ./hen intersection occurs, location 
of the M point at the intersection guarantees that both K v , bandwidth 
specifications have been satisfied, but may not guarantee the specifi- 
cation. In order that all three specifications be satisfied, the inter- 
section must lie in the area enclosed by the .3 curve and the B^ n axis 
in the first quadrant. 
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1.0 



Fig. 2-22a Solution of Example 2.5 using 

the 3rd order system band width cha,rt. 
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If one chooses the constant bandwidth curve for 6) bn =0.7 and 
draws the line with a slope of 0.35 from the origin as shown in Fig. 
2-22a, then the intersection point is at (.575, .204) and this point in- 
sures 5= *3 which is evident from the 173 curve configuration, thus 
one chooses the M point at M(a ln> a 0n) = MC • 575 , .204) to give 



t KfCt _ r~ n C” 

ln " Qo+KKif -■ S ^ 



a 0n “ 



K 



= . 24 # 



Since 6)^ : ,7 is chosen the normalizing factor a^ = 10-hKK a is de- 
termined from the relation 

=(/0 + K'ICa') 

from which 10hKK a = 28.6 is obtained. Substituting this value of the 
normalizing factor into the denominators of the two equations for the 
coordinate of the selected M point above, K,Kt and then K a are deter- 
mined and the solutions are 
K i 4675 
K t = 0.0987 

K a = 0.00388 
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2-6. Frequency response of the fourth order systems. 

The 17 curve of a fourth order system can be normalized so that a 
single curve can represent the 7 curve of any fourth order system. 
Suppose a fourth order characteristic equation is: 

F(S) = S 4 + 4* S%- 4, & -t 4 & = 0 

then the 7 curve is defined by 



(2-34a) 



| 3 | ~ cO 
Bo= 3.2 CO 2- - a 

Equation (2-34a) represents a concave down parabola whose axis of sym- 
metry is B^ = a. 2 ^/ 2 . with vertex at (a 2 a 3 / 2 , a^a^M). Bq : 0 occurs at 
CO r 0 and Ds a| . 

Assume a normalized frequency ci> tn such that 6)= ^tn t * ien 
Equation (2-34a) becomes 

B, = a a a 3 o> t * (2-35) 



2 

Let B ln = B 1 /a 2 a 3 and B Qn r B 0 /a 2 then Equation (2-35) can be writ- 
ten as 

[3 (r i = U3 tri (2-36) 

3oti - O)^- 0) t7 ? 

Equation (2-36) is exactly the same equation for the 7 curve as is asso- 
ciated with a characteristic equation whose coefficients of the S 4 , 

2 

and S terms are unity. 

If a curve defined by Equation (2-36) is constructed, then with 
proper scaling factors of coordinates and frequencies, it will represent 
exactly the same curve represented by Equation (2-34a) and thus 
curves of fourth order systems are normalized. 
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The normalizing property of the \ 0 curve for fourth order systems 
can be used to obtain the frequency response of any fourth order system. 



Suppose a fourth order system closed loop transfer function 



f($) s 



L‘ t & + 



(2-37) 



W K ^ i s“t-a 3 s 3 + + a,s + a* 

which has a second order polynomial in the numerator. For the case of a 
fourth order system with third order numerator, add b 3 S to the numerator 
of Equation (2-37). By using a transformation or frequency scaling 

S = if 4 

Equation (2-37) can be written as 






^2 2^ -f-b|5a .4 + a 0 






*0>=- 



ki&a. aS 



-J& -■ 

*■ t>i 3s. + ^3-2. 



■+ 3. 



(2-38) 



The systems represented by Equation (2-37) and (2-38) are equivalent 
since the transformation S = a|s is one to one. Note in Equation (2-38) 
that the coefficients of the and terms in the denominator are made 
unity, and also, the coefficients in the numerator and denominator of 
Equation (2-38), by comparing those of same power terms, are scaled by 

-3/2 

the same amount. For instance, b^ and a^ both are multiplied by a 2 , 
and likewise for a 2 and b 2 , etc. 

Let a 0n — 

a ln " <3* / 3-Z&3 
a 2n =a 2 /^. =i 

a 3n=5 3 /<2 3 =i (2-39) 

b ln = bt / <32% 
b 2n = \pz / & 2 - 
a 2t =a/V a 3 
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then Equation (2-38) is 



2. _ ' 1 " _ 

_C_/\_ bsn/l + Zzt (2-38a) 

+ a£*<s 3 + <?*>■» A + a-** 

By substituting s = j ^ , and noting that B ln r a),.^ B Qn 8 a) - (0 t „ 
which are the equations defining the normalized curve for fourth order 
systems, the magnitude of the frequency response is given by 






I ^r.7] 1 

| ~£ian^ ^ ^>tr> (&J7) j 



(2-40) 



thus the frequency response of fourth order systems is represented in 
terms of normalized quantities. 

The methods described in Sections 2-1, 2-2, 2-3 and 2-4 can be ap- 
plied to get the frequency response of the fourth order system. The 
frequency response of the original system represented by Equation (2-37) 
can be obtained from 






(dp-bzCO ) +j b,cO 

(a 0 -B 0 ) -4*j a; (a, - s,') 



(2-41) 



which results from applying Equation (2-6) to the system (2-37). The 
advantage of using normalized Equation (2-40) over the original Equation 
(2-41) is that, since the normalized curve for any fourth order sys- 
tem is a fixed curve, it is not required to plot the curve for each 
system. 

Knowing how to interpret Equation (2-41) on the Bq vs B^ plane which 
is associated with the original system, Equation (2-40) can be inter- 
preted on the normalized B^ vs B^ n plane without difficulty. 
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The quantities or coordinates in the normalized and unnormalized 
systems or Mitrovic's plane are related by Equation (2-39) plus the fre- 
quency relation 



1/2. 

a - 2 cO t7L 

The M point is located on the normalized B Q n vs B^ n plane by the normal- 
ized coordinate ( a ^ n » a Q n ). 

The imaginary term in the denominator of Equation (2-40) can be rep- 
resented in the same way as for Equation (2-41), namely by CO (a^ - B^) 
plot; however, it must be multiplied by a 2t to get the aj tn^ a ln~ B ln^ 
plot, and likewise for the imaginary term in the numerator of Equation 
(2-41) if a system of interest has b^ other than zero. 

In Fig. 2-23, the normalized V 0 curve for fourth order systems is 
plotted on the BQ n vs B^ n plane with the same scales for both axes. Any 
point Mn( a i n » a o n ) located on the Bq d vs B^ n plane is equivalent to the 
point M(apag) on the Bq vs B^ plane where a^ n , a^ and a^, a^ are re- 
lated by Equation (2-39). Any frequency C0^ n on the normalized Q curve 
is equal to CO divided by a^ where CO is the frequency on the 

unnormalized curve at the point which is equivalent to the 0) point 
on the normalized f 0 curve. 

The following table is a comparison of the V 0 curves for normalized 
and unnormalized systems; 



Axis of symmetry 

Coordinates of vertex 

Frequency at vertex 

Freq. at B^ intersect 

Coordinates of B^ 
intersect 



Normalized !7 curve 

Bm = 0-5 
(.5, .25) 

4) t n - .707 
"tn = o, 0) tn = 1 
( 0 , 0 ), ( 1 , 0 ) 



Unnormalized V 0 curve 
Bi = .5 a 2 a 3 

2 

( o 5 ^ 2^2 9 • 25 &2 ) 
COz .707 
CO = 0, CO = S a 2 
(0, 0), (0, a 2 a 3 ) 
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The equation defining constant bandwidth curves of fourth order sys- 
tems on the normalized Bq vs plane contains a parameter a2 ts , therefore 
it is not a practical way to try to provide a universal chart of constant 
bandwidth curves of fourth order systems , because for each a 2 t there re- 
sults a family of constant bandwidth curves. However, it would be worth 
while to investigate the effect of varying a 2 t on bandwidth, since if 
the effect of a 2 t on bandwidth is known, then if a constant bandwidth 
curve is drawn or bandwidth is determined for some &2t it may be 

possible to get some information about bandwidth for another a 2 t . 

Assume two fourth order systems such that when they are normalized 
they result in identical system equations except they have different 
a 2 t 's, namely a^. = k^ for one of the two systems and a^ t = k2 for the 
other where k^ and k2 are constants, k^4k2 . 

The two systems have a common M point on the normalized Bq vs 
plane. Suppose this M point is located as s.hown in Fig. 2-23. Since 
the normalized curve and M point are common to both systems 8 the 
6) tn (a^ n - B^ n ) plot will also be common, but each system will have a dif- 
ferent (a ln - B ln ) plot as shown in Fig. 2-23. (k^ = 1 and k 2 s 2 



are assumed for a simple illustration.) 

Since the bandwidth of a system is obtained at the point where some 
magnitude circle intersects the ~==? (a^ n - Bj^) plot, and noting that 

CLln 

the smaller a 2t causes the ^f=^ a ln ~ B ln) P lot to be displaced farther 
from the B ln = a^ n line than the larger a 2t does, any circle centered at 
the M point intersects the (a ln - B ln ) plot at lower frequency for 
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In its closed loop transfer function, then the magnitude circle to 
determine bandwidth has radius equal to J~2 - .0707 and it intersects 

with the two ^==(a^ n - B ln ) plots at P| and P 2 to give the bandwidths 



"bn = - 342 


for 


fat = 1 


"bn = -35 


for 


IK, -- 2 



It is obvious the above reasoning holds for systems whose closed 
loop transfer function has any number of zeros, thus a conclusion is ob- 
tained as follows; 

The effect of on bandwidth of a normalized fourth order system 
is to increase the normalized bandwidth with Increasing a 2t . 

It must be remembered that the normalized bandwidth is related with 

unnormalized bandwidth CO^ by 

As a simple application of the above conclusion, refer to Fig. 2-13. In 
Fig. 2-13, the hyperbolas, the constant bandwidth curves are for a system 

, , , 

This system can be regarded as some system which normalizes so that the 
coefficient of the S 3 term is unity with a^ t = 1. In other words, if for 
a normalized system 

C \ 

^“ + kS + S+ft* +2 - 

if constant curves are plotted with a^ t r 1» then there will result 

exactly the same configuration of curves as in Fig. 2-13. For the 
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M( a in > aQ n ) point as indicated on the curve C 5 , is known to be .5. 

If there is another normalized system with the same M point but with a 2 t 
slightly different from 1 , then C 0 ^ n of this system will not be far off 
from .5 in the well-known direction of deviation from .5. Same reasoning 
can be made if the M point is located anywhere on the vs plane,, 

The following numerical examples are to illustrate using the normal- 
ized 17 curve in obtaining frequency response of fourth order systems. 
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Example 2.6 



Given a closed loop system transfer function 

C fn\ / 6°0 

g* +22S 3 + I7ZS X + &90S + MOO 



(2-42) 



obtain the frequency response (magnitude). 

The problem can be worked out without normalizing the systems 5 how- 
ever, it will be shown on the normalized BQ n vs B| n plane. 

For the given system, the coefficients are 
a Q = 1600 
a L = 680 
a 2 = 172 
a 3 = 22 

By using relations between normalized and unnormalized quantities in 
Equation (2-39), the coefficients are normalized to 

s on = a 0 / aj 2 

a,* = & J & z 2l 3 = J8 

and fa 2t =fa x /a 3 W 

The frequencies of the normalized and unnormalized system are related 
by £0 =/aI lo) tn 



From Equation (2-40), the frequency response (magnitude) of the given 
system is, in terms of normalized quantities. 






■ oSV-Z 



In Fig. 2-24 the point M(a^ n , aQ^) = M *0542) is located and 

the lines B ln : a ln s .18, BQ n = a 0n = .0542, B 0n z (l + f~2 )a Qn s .131 
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are drawn on the normalized Bg vs B^ plane. 

Since the line B rt = (l+/l)a n = .131 crosses the Q curve at the 
Un On 

point where the frequency is slightly less than Ct) tQ - .4, and since the 
system has no zeros, the bandwidth 4)^ must be less than .4. From the 
V 0 curve and M point configuration on the normalized Bq vs Bjl plane as 
shown in Fig. 2-24, it is seen that the first one-half part of the 
curve is sufficient to see the frequency response of the system. 

In Fig. 2-25, the curve, M point and the B^ n s a^, Bg n r 
and Bg n = (1 + /T )ag n lines are drawn to expanded coordinate scales. 
Frequency response is represented in the same way as in Fig. 2-8, the 
only different thing from Fig. 2-8 is the i a ln “ B ln) pl°t» 



As illustrated previously, the (a-^ " B ln^ P^ ot con_ 

structed such that the horizontal distance from the line B lm = a 1r , s .18 
titn . 



D ln " rt ln 

to any point on the -r==r (a, - B. ) plot is equal to the horizontal 

distance from the line B^ n = a^ n = .18 to the curve multiplied by 



The — — (a. - Bi t ) plot bears the same frequency as 

/ft t 



tyn>. 

/ a 2t 

) 0 if the ordinates are the same as in the - B^) plot. 

Some magnitude circles are drawn to show the magnitude A 



f&zt 






< a ln" B ln) 



at some frequencies. The magnitude circle tangent to the 

(t)^ 

plot which is for resonant point touches to the ^ a ln“ B ln^ P^ ot: at 



C0 tn s .24 with radius equal to .05 to give 

the resonance frequency k) rn s .24 or 63 r = (13.1) 4) rn = 3.14 
the resonance peak = .0542/ .05 = 1.08 

The magnitude circle with radius equal to y^2 ag z .765 crosses the 

y=< a ln- E ln) P lot at “tn = • 



388 to give the bandwidth 
84 



W bn = .388 or 0) b s (13.1) (.388) = 5.09 
Substituting S = j&) in Equation (2-42), the magnitude for 6) s cd r r 
3.14 and Cd= (iC^ r 5.09 calculate to 

\%jW)\ = 1.1 for resonance peak Mp W 



and 



|^l : 



.71 for bandwidth 



Magnitude circles are drawn in Fig. 2-25 for several frequencies and mag- 
nitudes are indicates, thus the frequency response is obtained. 
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Example 2.7 

Given an open loop transfer function for a unity feedback control 
system as 



3A£(S+l){S+5} 



S (5+.333)(S+2 > )(5-H0 n ) 

determine the bandwidth of the closed loop system. 
The closed loop system transfer function is 
3/.6S z + t?t>S + 158 






(2-43) 



K' Jr s*+ i2.333S 3 +S$.6S*+/96.66£ + i-58 

which has a second order numerator. 

The quantities to be normalized to work on the normalized Bq vs 
plane are, using the relations in Equation (2-39), 
a Qn = 158/55.6 * .0512 
a ln = 196.66/(12.333) (55.6) = .286 
b ln = 190/(12.333) (55.6) = .277 
b 2n = 31.6/55.6 = .568 

= / 55.6/12. 333^ = .603 
■/a-) = 7 *44 
0) = 7.44 0) tn 

The frequency response of the system in terms of normalized quantities 
is defined by, using Equation (2-40), 






l-o 



(.05 -EU 



For the bandwidth 



+ j =/z\ (.m 2-.siS<hj (-^>1 
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The left side of Equation (2-44) can be represented by the usual 

( • 286 - B ln ) plot. On the right side of Equation (2-44), both real 
and imaginary terms are frequency dependent. As illustrated in Section 
2-3, the real term with the factor /”2 is represented on the B^ n = a^ n = 

.286 line and the imaginary term with the factor fl is represented on 

the BQ n r = .0512 line. Those are shown in Fig. 2-26. 

In Fig. 2-26, the normalized curve with the point M(.286, .0512) 
is shown. On B^ n = .286 line, /T b 2n^tn = *805 ft) is calibrated as 
follows. The point corresponding to C0 tn = 0 which is located below B^ n 
axis is at the distance of J~2 ag^ = .0723 from the M point. From this 
point up, the B^ n = .286 line is labeled with (3) tn such that the distance 

from the 0)^ = 0 point to a tt> tn = Ct) t ^ point is equal to 

2 2 

fl b ln = .805 4) tl so that at Q) tn r fl) tl , the distance from the 
M point to the 6) tn r £j) t ^ point on the B^ n = .286 line is equal to 

lyzCa^-b,^)! =|.^23-. serajj 

For example at Q) tn s .9 , J~2 (aQ n - b^ n 0) t jj) is equal to the dis- 

tance from the M point to the 6) tn s .9 point on the B^ n s .286 line, 
which is equal to .73. 

The imaginary term in the right side of Equation (2-44) is repre- 
sented on the BQ n s aQ n = .0512 line as follows. On the Bg n s aQ n s .0512 
line, each point to the right of the M point bears frequency a) tn such 
that a point corresponding to &) tn : on the Bg n = .0512 line is 

located at a distance of fl b^ n t j_/ J a£ t = .65&) t ^ . For example, 
at 6) tn : .9 on the Bg n = .0512 line, the distance from the M point is 
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(.65) (.9) = .585, thus the line B Qn = .0512 has linearly displaced o) tn 
whereas the line B^ n = .286 has quadratic c0 (;a . 

With those two lines, the quantity on the right side of Equation 
(2-44) can be read off at any frequency. For example, at = .9 the 

quantity /*2 j ( . 512 - .568 £0^) + J'g^C»277)j is the length of the diagonal 

of the rectangle formed by the point M, Od z .9 points on B, = .286 

tn J- 11 

line and BQ n = .0512 line, which is the radius of the magnituce circle 
C A2. 

Then the bandwidth is determined in the same way as in Example 2.3. 

In this example, as seen in Fig. 2-26, 

“bn -- ' 9 

Noting that frequencies are scaled by 6)= 7.44 &) tn , bandwidth of the 
original system is 

COb -- 6.7 

By substituting S = J 6.7 in Equation (2-43), the magnitude calculates 
to 

|-|(jb. 7 )|=.7b 7 
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2-7. Evaluation of phase angle. 

In the previous sections the magnitudes of the frequency response 
of a system are evaluated on the Bq vs plane. Since the numerator 
and the denominator of the equation defining frequency response can be 
represented as vectors on the Bq vs B^ plane, it is obvious that the 
phase angles in the frequency response of a system can also be evaluated 
on the Bq vs B-^ plane. 

In this section the method for evaluating the phase angles in the 
frequency response of a system is illustrated with some of the numerical 
examples which are considered in the previous sections of the magnitude 
evaluation. 

Again writing the general equation for the frequency response of a 
system as 

N( ]' 

^ ( a o- 8. ) -tj 

the phase of the frequency response is 



“H 




vS) = ^ M y ao-Bo) ,S) 



(2-45) 



Since the vectors N(jtO) and (a^ - B^) -*j'tO(a^ - B^) are measured 
from the M point and the real and imaginary parts of the vectors are 
represented on the Bq vs B^ plane in such orientations that they are 
parallel to the Bq and B^ axes respectively, a proper reference in meas- 
uring the angles must be set up. 

Consider the angle Z-((a Q - B Q ) -+j a3 (a 1 - Bq)) with the aid of 
Fig. 2-27. The A) (a^ - B^) plot is constructed as illustrated in the 
previous sections with the curve and the M point configuration as 
shown. Select a point such as P in Fig. 2-27 at frequency of GOj and 
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Fig. 2-27. Evaluation of the phase, angles 
in the frecpioncy response. 



line 



92 



draw the vector PM. The components of this vector define the real and 
imaginary terms (aQ - Bq) and ^(a^ - B^) in magnitude and sign. The 
horizontal component 6)j(a^ - B^) is the imaginary part and positive when 
directed to the right. The vertical component (a Q - Bq) is the real part 
and positive when directed upward. Considering the signs of the compon- 
ent vectors, the angle of the vector PM = (Sq - Bq) +j Cii^a^ - B^) is the 
angle 0^ which is measured from the vector (aQ - Bq) to PM, and this is 
equal to the angle measured from the line B^ = a^ below the M point to 
MP and positive when measured to the clockwise direction. The angle of 
the vector N(j^) is naturally which is measured from the line B^ » a^ 
above the M point to the vector ) with the positive direction as in- 

dicated, thus the phase angle of the system in the frequency response at 
frequency Ct)^ is 

It must be noted that the reference line of the angle of N(j£0 ) de- 
pends upon how the vector N(jCi) ) is represented and it is seen that the 
vector N(j 60 ) can be represented either to the right or left of the line 
B 1 = a l an< * e i t ^ ier above or below the line Bq = aQ ; however, for the 
angle of the vector (aQ - Bq) +j&) (a^ - B^), the reference line is always 
the line B^ = a^ below the M point and the positive angle is measured 
from the reference line to the vector MP in the clockwise direction, thus 
the 0°, 90°, 180° and 270° positions are as indicated. 

Evaluation of phase angles in the frequency response of an open loop 
system is also possible as the case of obtaining magnitude in the fre- 
quency response. The following examples illustrate evaluation of phase 
angles in the frequency response for closed loop systems as well as for 
open loop systems. 
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Example 2.8 



For the system equation in Example 2.1 which is rewritten below 




obtain the phase characteristics in the frequency responses for both 
closed loop and open loop systems. 

Closed loop system 

The closed loop system has the transfer function 



in which the numerator polynomial is of order zero and the angles 
/L N(J<W) = 0 for all frequencies, therefore only the angles 
^ ((Sq - Bq) +jw(a^ - B^)') are frequency dependent. 



from Fig. 2-8 which is for obtaining the magnitude of the closed loop 
frequency response of the same system. 

The reference angles 0°, 90°, 180° and 270° are indicated. The 
radial vectors directed from the M point to the <50 (a^ - B^) plot give 
the angles of the denominator of the closed loop frequency response equa- 
tion which is the vector (4 - Bq) 4-jCtf(2-B^) . Noting that the angles 
measured from the 0° reference line to the vectors are of the denominator 
and also noting that the angles of the numerator are zero degrees for all 
frequencies, the phase of the closed loop system at each frequency is the 
negative of the angle of the vector so measured, thus at Ci)= 0 the phase 
is zero degrees and as 0) increases the phase increases in the negative 
direction. For example, at 6) = 1.0 the angle of the vector measures 
45° and the phase of the system is -45°. The phase assumes -180° at 




In Fig. 2-28, the ^ curve and the ^(a^ - B^) plot are duplicated 
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Fig. 2-28 Evaluation of phase angles. 
(Example 2.8) 



4 B x 



CO = 1.4 and at this point the curve, ^(a^ - B^) plot and the line 
B 1 = a l = 2 intersect. If a protractor is overlapped on the figure such 
that the center point of the protractor coincides with the M point, then 
one can read off the phase at each Cd conveniently. 

Open loop system 

Rewriting the open loop transfer function as 

,S 3 + 35* + •zS 

one can construct the curve associated with the denominator polynomial 
of the open loop transfer function. By comparing the equation with the 
closed loop transfer function, it is obvious that the V Q curve for the 
open loop system is identical with that of the closed loop system. The 

Mq point defined analogously as the M point is, by noting that the co- 
o 

efficient of the S term in the denominator polynomial of the open loop 
system is zero, at B^ = 2 on the B^ axis and 

Mo=M 0 (s>o) 

is located in Fig. 2-28 as shown. Since the Mq and M points are on the 
same line B^ s 2, and the curves for the two transfer functions (open 
and closed loop) are identical, the <i)(a^ - B^) plot is not altered. 

By applying the same principles to the open loop system the phase 
can be evaluated with the vectors directed from the Mq point to the 

&) (a^ - B^) plot. Naturally the 90° 270° line must be shifted down to 

the axis. Noting that as GO tends to zero the ^(a^ - B^) plot ap- 
proaches the Mq point from the direction of the origin the phase at 
00 = 0 is -90° . As CO increases the phase increases in the negative 
direction. Some vectors representing the phase of the open loop system 
are drawn with dotted lines, for example at COr .5 the angle of the 
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vector measures 130° to give a -130° phase. 

At CO = 1.42 which is the point where the phase of the closed loop 
system is -180°, the open loop system also assumes -180° phase and it is 
interesting to note that this occurs for any system which has a frequency 
independent numerator polynomial. 

As 03 tends to infinity, from the behavior of the ^ curve and the 
CO ( a i - B^) plot it is seen that the phase approaches -270° which is the 
expected phenomenon. 

In the Bode diagram manipulation, the gain crossover frequency and 
the phase at that point are used for the stability check and as a figure 
of merit in a design problem. The gain crossover point which is the zero 
db or unity magnitude point can be found simply by drawing a magnitude 
circle with radius equal to Sq = 4 centered at the Mq point, which is to 

account the magnitude of the open loop system, thus the gain crossover 

o 

frequency is found to be 1.05 and the phase at this point is -168 to give 
+12° phase margin. This phase margin indicates a stable but too oscilla- 
tory system which checks with the curve and the M point configuration. 
The gain margin can be found by drawing another magnitude circle centered 
at the Mq point passing through the -180° phase point. 
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Example 2.9 

For the system in Example 2.2 

_ |a|P 00(5 + 47 ) 

J S (S I 6 70) 

the transformed closed loop system equation is, by using the transforma- 
tion S » 100s 

C,,N >*-!A ± 

In order to obtain the phase characteristic of the closed loop system 
the curve and the " B lt^ plot of Fis * 2-9 are re P roduced in 

Fig. 2-29. 



The angles of the denominator vectors 0^ are measured in a similar 
way as in the previous example. Noting that the system has the first 
order numerator and the vector N( j&> t ) at each 6) t is represented by a 
vector directed from the M point to the line Bg t = 2a^ t (the two vectors 
for 0> t = .2 and C0 t = .3 are shown in the figure with dotted lines) the 
angle Z.N(j’Ci) t ) at each is measured from the line B^ t = a^ t above the 
M point to the vector N(j’£D t ) at each b) t . For example at 
C0 t = ,2 } Zn(j .2) = 29°. Denoting the angles of the denominator and 
the numerator vectors at the same frequency by 0^ and 9^ respectively, 
then the phase of the closed loop system s “ ^1 » ^ or 

example at 6) t = .2, Z.A(j.2) = 29° - 31.5° - -2.5°. 

For the open loop system the transformed system equation is 

13.1 4 + 4.9-2 

A* + .t6 14 

The 1 curve for the open loop system is not changed from that of 
the closed loop system but the Mq point is moved to the point (.167, 0) 
which is seen from the coefficients of the s^ and s® terms in the 
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denominator of the open loop transfer function. This implies that a new 
G) t (a^ t - b k-) plot for the open loop system must be constructed with 
a^ t = .167 in order to obtain the frequency response of the open loop 
system and the procedure to construct this plot is straightforward. 

Since the Mq point is so close to the origin, the <^(.167 - B^ t ) plot 
will not give good resolutions of the angles of the denominator for 
small COj. if the plot is constructed in Fig. 2-29. This can be overcome 
by plotting the curves for small o) t on an expanded scale. 
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3. The fourth order system charts. 

In references 1 and 2, a universal third order system chart is con- 
structed on the Bq vs plane so that a single chart will provide solu- 
tions to any third order system problem. 

This is done by using a suitable transformation. Any third order 

system is transformed into a normalized system so that it has unity co- 
3 2 

efficients of the terms S and S in the characteristic equation to give 
a set of fixed Bq vs B^ curves. 

In a fourth order system, assuming unity as coefficient of the 

tern in its characteristic equation, Bq vs B^ curves are changed by chang- 

3 

ing either a^ the coefficient of the S term or a£ the coefficient of the 
2 

S term. 

Since for any polynomial it is possible to make any two of the coef- 
ficients unity (or the same numbers), by a suitable linear transformation, 

in a fourth order system characteristic equation any two of the coeffi- 
4 3 2 

cients of S , S and S terms which are the coefficients affecting the 
Bq vs B^ curves can be made unity so that all Bq vs B^ curves of any 
fourth order system depend upon only one coefficient. 

For each value of the variable coefficient, there results a different 
set of the Mitrovic's curves. If the sets of the curves for some suitable 
pertinent values of the variable coefficient are constructed and if the 
behavior of the set of the curves for the variable coefficient changing 
from one pertinent value to others is understood, then the sets of the 
curves may serve as fourth order system charts which can be applied to 
any fourth order systems. 



101 



3-1. Transformation of the fourth order system equation. 



Suppose a fourth order system characteristic equation is: 

F (s) = + a,;s% a, s 1 •+ <3_ 0 

and associated Bq vs B^ curves are defined by 

B 1 = 

B 0 -- 

Using a transformation or frequency scaling 

S=a 3 ^s 

Equation (3-1) can be written as 

F ( <f?3 A) = * "+■ ^3^ ■+■ % ^3 ^ "+■ 2, -4- a o — o 

or by dividing through by Q. .j 

£(/0 = /!*+ ^ 3 + 3./*/ 4- 

or (/$) — /6 ^ + yd 3 ■+■ <3^^ /S + 2. a t = 0 (3-3) 

where = -2. 2 / aj 2- 

a (t =c2, / 2j 

a 6t =a B / a/ 

Thus Equation (3-1) is transformed into Equation (3-3) in which the 
two coefficients of the and S'* terms are unity. 

Equations (3-1) and (3-3) are equivalent since the transformation 
S = a^s is linear and one-to-one. 

The Bq vs B^ curves associated with Equation (3-3) are defined by 

B| t = Z2t&(-S )&)„r + <h tint 

*. (3-4) 



(3-1) 



(3-2) 
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Equations (3-2) and (3-4) are equivalent to each other with the fol- 



lowing relation* 

=4)^/a 3 

Comparing Equations (3-2) and (3-4) which are Mitrovic's equations 
for a system represented on different coordinate scales, it is seen that 
studying a fourth order system with the second equation is more convenient 
than with the first, because in the second, only one variable coefficient 
ajj. appears whereas in the first the two coefficients a 2 and aj appear. 

A fourth order system closed loop transfer function represented as 



-)- bi£ 

a,£ + a* 



(3-5) 



transforms into 



C, .V -t- Act 



(3-6) 



to which the Mitrovic's curves defined by Equation (3-4) associate. 

The coefficients and frequencies of the transformed system are related 
with those of the original system by 
k/&3 
kt~ k/2j 
%= 

plus the relations defined as in Elation (3-3) for the coefficients of 
the denominators. 

The velocity constants, Ky and of the original and transformed 
systems respectively are related by 
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fee -I <v/*3 



since 




Thus a fourth order system is transformed into a system from which 
resulting Bq vs curves are dependent upon only one coefficient a 2 t . 
Comparing Equations (3-6) and (2-38a), it is noticed that the characteris- 
tic equations are of different forms which indicates that the normalized 



Equation (3-4) and any I ^ curve defined by Equation (3-4) cannot be 
plotted on the normalized Bq vs plane unless the system Equation (3-6) 
is transformed into the form of Equation (2-38a) and Equation (3-4) is 



can be plotted on the normalized Bq vs B^ plane; however, appears 

in too many terms and more complicated behavior of the curves is expected. 
Thus the normalized Bq vs plane is not suitable to study fourth order 
systems. However the normalized curve can do an excellent job in determin- 
ing absolute stability of a fourth order system. 

The coefficient a2t appearing in Equation (3-4) indicates that Bq vs 
B^ curves for a fourth order system depend upon the value of a 2 t , 



Vo curve in Section 2-6 is different from the curve resulting from 



forced into a form to give the normalized V curve. If one wishes to work 
on the normalized Bq vs B^ plane with ^ curves, then the system Equa- 
tion (3-5) or (3-6) is transformed into the form olT Equation (2-38a) to 



give the equation defining Mitrovic's curves as 



B ln * f&it -+• 5 

»0n = a '*‘ 



(3-7) 



Equation (3-7) gives the normalized Vq curve, and the 



curves 
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therefore it seems worth while to investigate the possible range of &2t° 

The coefficients a 2 and a 3 of Equation (3-1) » the characteristic 
equation of the original system, can be written as 

a 3 s r i* r * + r 3 + r * 

(3-8) 

a 2 = r , r * + r . r 3 + r ) r * + r a r 3 + + r 3 r * 

where -r^ , -r 2 , -r 3 and -r^ are the roots of Equation (3-1). 
The roots r^ , r 2 , r 3 and r^ fit into one of the following three cases; 

1) All 4 roots are real. 

2) Two real and one pair of complex conjugate roots. 

3) Two pairs of complex conjugate roots. 

2 

Noting that a ^ - a 2 /a^ , by using a 2 and a 3 expressed in terms of 

the roots, it is possible to determine the range of a^ fc to cover most of 

the fourth order systems for Case 1 and Case 2. 

Case 1 

„ _ r, r* + r, r 2 +■ n r u + r* r 3 + r 2 r* + r 3 r* 

2t “ (r, + r * + r 3 4-u) a 

nn r.r 3 t-rrr»tr a ifs-»-r a r» 4- 

rr+r/ +r*+rf + *( r,r n r 3 + r iv + r 3 $+ & r* + r 3 r* ) 

noting that r^+r/t r/ + >0 

o < <.5 

It can also be seen that for a fourth order system, if its open loop 
transfer function has a numerator polynomial of order less than two, and 
if it has all real poles, then 

a zt <.5- 

In such a system, the coefficients a 2 and a 3 of its characteristic 
equation can be expressed in terms of the poles of the open loop transfer 
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function. Those poles of the open loop transfer function replace the 
roots of the characteristic equation in Equation (3-8) to get the same 
result. 

Case 2 

Let the two real roots be r^ and ^ 9 and the complex conjugate 

roots r. and r, . 

3 4 

r 3 = 

r 4 = 5 - j 

and = ( r, -+ r a + a. -s +njV 4 5^ fr,+rO -+• 43X? 

a 2 = r,r 2 + 

to compare and a 2 

r,+r l r i +r i %25^(r,+r^4(4/-i)< (3-9) 

Assuming a stable system, Equation (3-9) is examined for several 
situations. 

If r l > r Z > ^ then 

3^- 3^ > (ll 5*- implies that 

if % >• 3°^. then > 0 and 3^ < \.Q 

It must be noted that the inequality above may have considerable tol- 
erance which implies that even for % some value less than .302 ag,. may 
be less than 1. 

Let r£ > r^ and by neglecting r^ in Equation (3-9) 

a^-a A > ^-+-2^60,,^ + (4 3 2 - \)0X? (3-10) 

For a well behaving system, t 2 is usually greater than the 

real part of the complex conjugate roots. 
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Let T 2 = k^&O^ where k is a positive real number greater than 1. 
Substituting this r 2 into inequality relation (3-10) 

^-2* >((£*+2^ iO^ (3-11) 

Inequality relation (3-11) indicates that for a reasonable number k, if 
£ is not too small, then a| - a 2 is greater than 0 and a 2t <1.0. 

For example if k=2.5 then a| - a^> 0 for £>.25 and if k-=3 
then a| - a 2 > 0 for t > . 229 

Thus for a fourth order system, if it has no more than one pair of 
complex conjugate roots of the characteristic equation and if it is ade- 
quately damped, then the coefficient a^ t is less than 1. 

For Case 3, which is the case of two pairs of complex conjugate roots 
present in the characteristic equation, it is hard to set up an upper 
bound of a 2 t as in Case 1 and 2. This can be seen as follows s 
For Case 3 

* rk ) = fV 

n = 5 A, -j 1 G = % Olna -j ^ 

(2 5, 0) 7i i "+ -2 % 

^2, 58 C&n i -+• ^ •+" ^a. 

<3 a . 

^711 ^ ^2 <2- 7 ^ 

Jb_ I + +#* 

a l ~ * 5* + <?5, Sa £ + 



Let 

and 

then 

and 

Let 

then 
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' 



Thus the ratio a 2/ 2 
a 3 

and $ 2 are sma H» then a 2/ 



depends upon ^ ^ 



2 assumes a large number. 
a 3 



k . If both $ i 
However, if a well 



damped system is assumed, then Case 3 can be approximated to Case 2. and 

a 2 t < 1 is a reasonable upper bound. For example if = .3, $2 s .7 

and $L z (jd 7)X then = .71 

a2 4.0 



In this section, in studying fourth order system curves, a£ t < 1 is 
assumed. It must be noted that a 2 t <1 is not the sufficient condition 
for a fourth order system to be well damped. Even though a 2 t is very 
small, a fourth order system may be unstable depending upon the location 
of the M point. It turns out that one can find a necessary condition on 
a 2 t for a system to have a relative stability greater than a given 

system damping coefficient. This will be discussed later. 
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3-2. Three dimensional representation of fourth order curves. 

The Bq vs B l curves for fourth order systems defined by Equation (3=4) 

2 

depend upon a 2 t which is the coefficient of the s term in the character- 
istic equation of the transformed system equation. a^ t can assume any 
value within its allowable range and independent of frequency on the Bq vs 
B^ plane. 

If a new variable B 2t = is defined and B 2t varies on B 2t axis 
which is chosen to be perpendicular to both B^ t and Bq c axes to form a 
three-dimensional cartesian coordinate system, the equations 

Bzt = 

B, t = + ^ (3-12) 

B. t = a Jt 

define the Mitrovic's surfaces on the B 2 t - B^ t - Bqj. coordinate system. 

For fourth order systems. Equation (3-12) defines fixed surfaces 
which can be applied for any fourth order system as the case of the third 
order chart for third order systems. 

Since the variable B 2 t is continuous and Bq c vs curves are con- 
tinuous, the Mitrovic's surfaces are also continuous. Noting that any 
plane B 2 t = s constant which cuts the surfaces gives the set of BQ t vs 
B lt curves for various values of ^ , by applying the continuous property 
of the Mitrovic's surfaces, all of the basic property of the Bq vs B| 
curves and the analysis and design principles of feedback control systems 
on the Bq vs B^ plane can be extended to the three-dimensional Mitrovic's 
surfaces. 

For example, the stable region for the M point on the Bq vs B^ plane 
which is the region in the first quadrant enclosed by the curve can be 
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extended to the three-dimensional Mitrovic's surface as the region (volume) 
in the first octant enclosed by the surface which is the surface made 

up by the curves for all values of > 0. Locating the point 

M(a lt> a 0t ) on B Qt vs B lt plane is to locate the point M(a 2 t , a lfcS a 0t ) in 
the Mitrovic's volume. 

Fig. 3-1 is a schematic sketch of the surface. B 2 t axis is chosen 
so that B lt an< * B 0t 3X68 f Qrm a right-handed coordinate system. The 

equations defining the surface are 

(3=13) 

Sot = " ^71 1 

which are the equations for a parabolic surface. The surface forms a 
cave with closed end at the origin on the Bq^. * 0 plane. Applying the 
stability criterion to each BQ t vs B^t curve which is formed by cutting 
the surface with each B 2 t = a 2 t plane, the stable region is the whole vol- 
ume inside the cave. 

If the surface is cut by a plane 3 q c = aQ t = constant, then there re- 
sults a B^ t vs B£ t (or B 2 t vs B^) curve for £ s 0 on the plane and like- 
wise for the Bqj. vs B 2 t curve. 

It is interesting to note that, Bg t vs B 2 t curves for $ s 0 on Bq^. 
vs B 2 t planes for fourth order systems are straight lines whose slopes and 
axis intersections are determined by a^ t only and an extremely simple 
method to test the absolute stability of a fourth order system can be 
derived. 

By eliminating fiO nt and a 2 t in Equation (3-13), the surface of a 
fourth order system is 
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Fig. 3-1 The surface for 4th order systems. 
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By cutting the surface represented by Equation (3-14) with a 
B lt = a lt Pl ane > the B 0t vs B 2t curve for ^ ^ is obtained as 

B ot = (3-15) 

which is an equation defining a family of straight lines on the Bg t vs 
B 2 t plane (or B^ t z a^ t = constant plane). 

For each value of a^ t , there results a different line. By reason- 
ing with the stable region in the three-dimensional volume, the region 
for the point M (a 2 t , ag t ) on Bg t vs B 2 t plane to guarantee a stable 
system is the area in the first quadrant below the line defined by Equa- 
tion (3-15) or mathematically 

0< a„ < a.t • a* t -a, t 

which agrees with the Routh criterion. 

The r 0 curves for some values of a^ t are plotted in Fig. 3-2. As 
the Equation (3-15) indicates, a line for a^ t = k cuts the B£t axis 
at B 2 t = k and it has a slope of k. The frequency on these lines is 
constant at all points on a line which is a peculiar nature of the 
curve on a Bg t vs B 2 t plane. Thus the line itself is a constant k) nt 
curve on the Bgt vs B 2t Pl ane * 

For instance, along the line for a^ t s .2, = .446 which is 

equal to as the second of Equation (3-13) implies. Thus a point 

M (ag t , a 2 t ) located anywhere on a single line, say for a^ t = a^ fc , 
represents a closed loop system which has a pair of pure imaginary roots 
r = = J V a it which is readily seen from the basic nature of the 

curve. 

As an example, suppose a characteristic equation of a fourth order system 

—/6 U +/4 2 ■+• +■ *2"<J -+ 
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Since a^ t = .2 the line is the line labeled as a^ t = . 2 a 
CO t = .446 in Fig. 3-2. The Mitrovic's working point M (aQt » a 2t) 
cated anywhere above the ^ line represents an unstable system^ and 
dually for below the line. If the M point is followed along the line 
then resulting every characteristic equation has a pair of imaginary 
roots r s ±j .446. For instance, if M (.7, .1) is chosen, then the 
characteristic equation is 

U)= A** A 3 +-'7A*+ .2/ + .l=0 

and -f,(±j.w6)-=0 

The point M (.7, .1) is also on the line for a^ t : .5, therefore 
another characteristic equation 

= +. 7 / 1 * +. 5 * + .|=0 

has the roots r r ±j/lf and it calculates to 

£.$/?)= o 

From the Mitrovic's basic equation and the way of defining a three- 
dimensional plot, the concept of a Mitrovic's surface can be generalized 
to any order system and to a set of any three coefficients of the charac- 
teristic equation. Then fixing one of the three axes to a constant, 
there results a well-defined curve which is the one defined by the gener- 
alized Mitrovic's equation in Reference 3, though the generalized Mitro- 
vic's equation turns out to be undefined and uninterpretable at some 
particular values of ^ for some particular set of coefficients. (For 
example, the curve for the BQ t vs B2t plot.) 

The most significant feature in the concept of the Mitrovic's sur- 
face is the possibility of controlling three coefficients simultaneously 
since the surface is independent of the three coefficients which are to 
be controlled. 
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3-3. Study of fourth order curves. 



As pointed out in Section 3-1, a fourth order system can be trans- 
formed into an equivalent system so that there remains only one param- 
eter (coefficient) in the characteristic equation on which the fourth 
order system curves depend. This feature greatly simplifies studying 
the behavior of the Mitrovic's curves for fourth order systems. 

In Section 3-2, it was shown that if a fourth order system charac- 
teristic equation is transformed such that the coefficients of the 
and S'* terms are unity, then the coefficient of the S^ term a 2 t which is 
the only one affecting BQ t vs B^ t curves is less than 1.0 for most prac- 
tical fourth order systems. 

Since the BQ t vs B^ t curves are continuous for all values of ^ nt s 
^ , and the coefficients of characteristic equation, if the behavior of 
the curves is well investigated and if the curves for some pertinent 
values of parameters are constructed, then those preconstructed curves 
may serve as a guide or charts for both analysis and design of fourth 
order systems 

In this section, the effect of a 2 t on Bqj. vs B^ t curves is studied 
to provide the fourth order charts. 

In Equation (3-4) which is the equation defining Bg t vs B^ t curves 
for a fourth order system, there are two parameters a 2 t and $ . For one 
value of a 2 t , each gives a different curve to form a family of 
curves corresponding to the particular value of &2t’ Th *- S family is Re- 



sults a different family. 

Fig. 3-3 to Fig. 3-22 show such families. On each sheet the curves 
for five values of £ are shown. For one value of a2t» there are 



noted by 




or simply 




Then, for each a 2 t there re- 



115 




curves(4th order systems) 
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Fig. 2-5 B vs B , curves (4th order systems) 
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Fig. 3-14 B Qt vs B.^ curves(4th order systems) 
for a =.6 (5=.5, .6, .7, .8, .9,1.0) 
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I’ig* 3-18 B Qt vs B lt curvos(4th order ayntema) 
for a =.8 (S=.5,.6,.7,.8 f .9,1.0) 
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associated two sheets of graphs, one for 0 ^ "5 ^ • 5 and the other for 
.5 6 5 ^ 1.0 with interval of .1. The values of a^t are from .1 to 
1.0, thus for every two sheets of graphs, a 2 t increases by interval of 
.1. (Note that for .4 ^ a 2t~ t * ie fi rst graph for each value of a 2 t 

contains the family for 0 ^5^ 1.0 and the second graph con- 
tains the family for -5 £ £ i 1.0.) The constant 0) nt curves 

are also shown for some pertinent values of &) nt * 

It is observed that as a 2 t increases , the curves are in general ex- 
panded in both abscissa and orxinate. This is well expected from the 
equations defining Bg t vs B^ t curves. The equations are rewritten below 

B lt + 

, u. (3-4) 

3 0t =a 2 

The terms with which a 2 t contributes to B^ t and BQ t are non-negative 
for 0 — 1.0, thus increasing a 2 t increases Bj^ t and Bg t . 

For small values of a 2 t each member of the family is so 

shaped that the whole region under the curve in the first quadrant is 
the region for the M point where the relative stability is guaranteed to 
be greater than the value of £ of the curve. Such a region is denoted 
by R^ . For example in Fig. 3-3, the whole area enclosed by the ^ 

curve and B^ t axis in the first quadrant is the region R # 5 on which if 
the point M (a^ fc , a g t ) is located, then the system has greater than .5. 

As a 2 t increases, the manner in which the R^ region is formed is 
changed in such a way that the curves for high ^ do not form R^ . 

The behavior of the curves depends upon the values of a 2 t and ^ . 
Examination of Equation (3-4) reveals that for a particular value of 
a 2 t , there is a certain critical value of ^ such that if is increased 
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beyond the critical value, then the curve exhibits counter-clockwise 



near to the origin or no R^ at all. Equivalently for a particular 
value of ^ there exists a critical value of a£f 

Investigation of the critical value of ^ or a 2 t is important since 
in synthesis of a system one must choose a point M (a^ t , ag t ) so that the 
specified system performances are met and particularly for the relative 
stability of the system, the M point must be located on a R^ where ^ 
is some suitable value to meet the specified relative stability. For 
example if the system relative stability is specified by ^ = .5 then 
one should not choose the M point anywhere on the curve for 

a 2t = 1.0. (See Fig. 3-21 and 3-22.) 

In synthesis of a system, usually is chosen at some pertinent 
value according to the required relative stability. Therefore it is con- 
venient to investigate the critical value of a 2t for given values of 

, thus a study of the behavior of the curves with changing a 2t for a 
fixed $ is needed. This requires that the family of curves must be 
represented in a different fashion from those of Fig. 3-3 through Fig. 
3-22 in which the families of Bqj. vs B^ t curves with changing ^ f° r 
fixed a^ t (the family (^^^ a 2t) are re P resente ^* That is s the family 

of curves must be represented with changing a 2( . for fixed ^ . Those 



At this point, in connection with the study of the families 
1 £){ a 2t\ » a constant €> plane is defined as follows; 

A constant ^ plane is a Mitrovic's plane on which Mitrovic's curves 
for fixed 5 are plotted, and denoted by the ^ = k plane where k is a 



encirclement on the Bg t vs B^ t plane so that there is formed a small R ^ 



families are denoted by 
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